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Preface 

\, 

'A  model  of  a  system  is  a  simplified  version  of  the  real 
system.  The  simplification  is  done  by  introducing  a  set  of 
assumptions  that  express  our  understanding  of  the  system's 
behavior.  In  the  process  of  simplification,  we  omit  non-dominant 
effects,  leaving  the  main  features  of  the  system.  The  objective 
of  simplification  is  to  obtain  a  set  that  can  be  solved  by 
available  tools  and  still  maintain  the  main  charact erist ics  of 
the  real  system. 

The  model  is  a  tool  that  provides  forecasts  of  the  response 
of  the  real  system  to  various  excitations.  Thus  it  is  essential 
that  the  model  retains  those  features  that  are  the  subject  of 
invest igat ions . 

This  volume  contains  7  papers  dealing  with  models  of  brain 
mechanics.  The  objective  of  the  models  is  to  provide  a  tool  for 
simulating  the  mechanical  behavior  of  the  cerebral  system,  as 
manifested  by  pressure,  velocity,  stress  and  strain  variations, 
in  response  to  changes  in  input  pressures  and  fluxes.  ■(- - 

At  this  stage,  the  models  deal  only  with  changes  in  fluxes 
pressures,  stresses  and  strains.  At  a  later  stage  of 
investigations,  the  transport,  spreading  and  interactions  with 
the  tissue  material  of  chemicals  carried  with  the  brain  fluids, 
will  be  investigated  and  reported. 

In  principle,  our  objective  is  to  develop  continuum  models 
of  brain  mechanics.  Such  models  will  provide  information  on  the 
behavior  at  every  point  in  space  and  time  of  the  brain  domain,  in 
response  to  changes  in  the  normal  behavior  in  various  parts  of 
the  brain  (e.g.,  vessel  occlusions).  However,  at  this  stage 
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reported  here,  only  one  simplified  model  that  investigates 
deformation  and  stress  in  response  to  pressure  variations,  is  of 
the  continuum  type.  All  other  models  are  of  the  multicompartment 
type . 

In  a  mult icompart ment  model,  the  entire  domain,  here  the 
portion  of  space  occupied  by  the  brain,  is  visualized  as 
comprised  of  a  number  of  interacting  compartments.  Each 
compartment  reprsents  the  lumped  behavior  of  a  certain  part  of 
the  brain  that  has  distinct  features.  Fluids  move  from  one 
compartment  to  the  next  under,  pressure  gradients.  In  the 
present  investigations,  a  seven-compartment  model  was  employed  to 
simulate  the  behavior  of  the  various  parts  of  the  brain. 

Every  model  contains  a  number  of  parameters  that  express  the 
excitation  response  characteristics  of  elements  of  the 
invest  igated  system.  In  order  to  make  use  of  a  model  in  specific 
cases,  the  values  of  these  parameters  must  be  known.  The  only 
way  to  determine  these  parameters  is  to  make  use  of  observations 
in  the  real  cerebral  system,  and  compare  them  with  values 
predicted  by  the  model.  The  process  of  determining  model 
parameters  is  referred  to  as  the  inverse  problem. 

Accordingly,  one  of  the  papers  describes  how  the  basic 
parameters  of  the  model  are  estimated. 

f 

Typically,  pressures  and  fluxes  in  the  brain  vary  in  a 
cyclic  manner,  occasionally,  especially  under  pathologic 
conditions,  transient  changes  may  take  place.  Accordingly,  the 
models  presented  in  the  various  papers,  simulate  non-steady  and 
transient  flow  conditions.  In  order  to  enable  such  conditions, 
compliances  had  to  be  introduced  between  adjacent  compartments 
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(in  addition  to  resistances  that  control  the  pressure  flux 
relat ions) . 

Results  of  model  predictions  were  compared  with  actual 
observed  data,  whenever  available. 

The  results  of  the  comparisons  were  satisfactory  and 


encouraging . 
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INTRODUCTION 


Ever  since  the  cerebrovascular  fluid  system  was  modelled,  the  rela¬ 
tively  simple  compartmental  representation  has  been  favoured.  MONRO,  the 
younger  (1783),  the  pioneer  of  intracranial  mechanics,  characterized  the 
physical  forces  in  the  intracranial  cavity  assuming  incompressibility  of  the 
fluid;  nearly  incompressibility  of  the  brain  tissue,  and  zero  motion  of  the 
boundaries.  The  intradural  space  was  regarded  as  bi -compartmental  -  brain 
and  blood,  so  that  a  change  in  either  compartment  had  to  be  compensated  for 
by  the  other. 

KELLIE  (1824)  modified  this  hypothesis  by  assuming  three  instead  of  two 
compartments,  namely  arteries,  veins  and  brain  tissue.  About  the  brain  tissue, 
Kellie  said:  "The  brain  itself,  little  compressible  is  contained  within  a 
firm  and  unyielding  case  of  bone  which  it  exactly  fills  ....  If  these  premises 
be  true,  it  does  not  then  appear  very  conceivable  how  any  portion  of  the  circu¬ 
lating  fluid  can  ever  be  withdrawn  from  within  the  cranium,  without  its  place 
being  simultaneously  occupied  by  some  equivalent  or  how  anything  new  or  exu¬ 
berant  can  be  intruded,  without  an  equivalent  displacement"  (loo.  cit.,  p.  102). 

The  Monro-Kellie  doctrine  of  almost  absolute  rigidity,  prevailed  into 
this  century  and  was  only  relaxed  in  stages.  The  number  of  fluid  compart¬ 
ments  was  increased  to  six:  artery  (A),  capillary  (C) ,  venous  IV),  venous 
sinus  (S) ,  jugular  bulb  (J)  and  cerebrospinal  fluid  (F)  (AGARWAL,  1971).  Yet 
the  fluid  itself  remained  incompressible.  More  recent  approaches  relaxed  the 
latter  assumption  and  the  fluid  was  taken  as  linearly  compressible,  namely 
the  change  in  pressure  and  the  change  in  volume  are  proportional  to  each 
other,  and  the  coefficient  of  proportionality,  the  bulk  modulus  or  its  inverse 
-  the  compliance,  are  constants. 
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The  thick  layer  of  the  brain  tissue  between  the  ventricles  and  the  sub¬ 
arachnoidal  space  has  also  undergone  various  stages  of  modelling.  The  buoy¬ 
ancy  theory  (LIVINGSTON,  1965)  regarded  the  tissue  phase  submerged  in  the  ideal, 
incompressible  (Pascal ian)  CSF  fluid  to  which  Archimedes  law  applies.  The 
tissue  itself  was  taken  as  single-phase  and  incompressible. 

When  the  incompressibility  assumption  was  abandoned,  it  paved  way  to 
compressible  representations  for  the  brain  tissue.  Numerous  experimental 
studies  have  been  conducted  in  the  past  four  decades  in  that  direction.  The 
results  pointed  at  the  inelasticity  of  the  tissue  or,  alternatively,  the  non¬ 
linearity  of  the  compliance  (RYDER  et  al.,  1953;  MILLER  and  GARIBI ,  1972; 
LUNDBERG  et  al.,  1974;  MARMAROU  et  al.,  1975;  MILLER,  1975;  HAKIM  et  al.,  1976; 
LEWER  ALLEN  and  BUNT,  1977;  GRIFFITH  et  al.,  1978;  BRUCE,  1978;  CHOPP  and 
PORTNOY,  1980).  To  overcome  the  "non-linearity"  of  a  single  coefficient,  more 
complex  models  of  single-phased,  multi-parameter  viscoelastic  materials  were 
introduced  such  as  the  one  consisting  of  four  viscoelastic  coefficients  known 
as  the  "three-parameter  solid",  coupled  dynamically  with  another  elastic 
element  (PAMIDI  and  ADVANI,  1978). 

In  the  literature  on  the  mechanics  of  head  impacts,  the  skull  was  ide¬ 
alized  to  be  a  rigid  sphere  with  an  opening  that  simulated  the  foramen  magnum, 
and  the  spinal  dura  mater  was  idealized  as  a  cylindrical  membrane  fitted  to 
the  foramen  magnum  (LIU,  1978).  The  intradural  content  of  the  central  nervous 
system  (CNS)  fluid-filled  continuum  was  regarded  by  some  as  a  single-phase 
incompressible  fluid  in  a  two-compartment  (skull  and  spinal  dura  mater)  struc¬ 
ture;  or  else  as  a  single-phase  compressible  elastic  fluid  (later,  also  a 
viscoelastic  fluid)  possessing  a  single  averaged  "elastance"  (strictly,  a  bulk 
modulus  -  inverse  of  compliance)  along  with  a  shear  modulus  (POLLACK  and 
BOSHES,  1936;  OMMAYA  and  HIRSCH,  1971;  LaFGREN  and  ZWETNOW,  1973;  GOLDSMITH, 
1972;  MARMAROU,  1973;  KING  and  CHOU,  1976). 


Here,  we  extend  the  discussion  to  the  n-compartmental  model  and  present 
the  methods  how  to  solve  it  for  the  general  non-steady  state  flow  with  con¬ 
stant  resistances  and  compliances.  An  explicit  numerical  solution  is  given 
for  the  case  n= 7.  It  is  also  shown  that  for  the  slow  mode,  compartmental  pres 
sure  waves,  solutions  for  the  more  simplified  quasi-steady  state  flow  are  al¬ 
ready  in  good  agreement  with  some  clinical  measurements. 

1.  The  Compartmental  Flow  Equation 

By  its  nature,  the  compartmental  modelling  of  any  hydraulic  system  is  a 
lumped-parameter  modelling.  Thus,  the  resistance  to  flow  due  to  a  particular 
vessel  type  is  lumped  at  the  outflow  of  the  compartment.  Likewise,  the  inte¬ 
grated  change  in  volume  of  each  compartment  is  representable  as  an  overall  com 
partmental  property.  Furthermore,  if  there  is  a  production  or  a  drainage  of 
fluid  in  a  compartment,  the  source-sink  function  is  attributed  to  the  entire 
compartment.  Finally,  the  functional  interaction  between  the  components  of  a 
lumped  parameter  model  is  assumed  to  be  at  the  interface  between  adjacent  com¬ 
partments. 

The  general  compartmental  flow  problem  is  that  of  the  non-steady  state. 
Here,  all  parameters  are  functions  of  time  and  the  deformabi 1 i ty  of  the  com¬ 
partments  is  taken  into  account.  Fluxes  are  composed  of  two  terms:  (i)  Flow 
rate  of  fluid  as  result  of  pressure  difference  (transmural)  between  two  com¬ 
partments,  expressible  as 

AP 

q>  =  ZAP  =  —  (P  =  pressure,  ft  =  resistance,  (1) 

n 

z  =  fluidance) 

(ii)  Flow  rate  from  the  deformational  volume  change  which,  in  turn,  relates  to 
the  change  in  pressure  by  a  functional  relation  of  the  type 
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Q"  =  =  C  (C  =  compliance)  (2) 

The  governing  equations  for  the  compartmental  flow  problem  are  the 
balances  of  mass  and  of  linear  momentum.  We  do  not  consider  high-speed  turbu¬ 
lent  motions  (vanishing  of  the  angular  momentum)  and  the  system  is  assumed 
isothermal  and  thermodynamically  stable  (identical  vanishing  of  the  energy 
balance).  Under  these  assumptions,  the  mass  and  linear-momentum  balance  con¬ 
ditions  for  the  n-compartmental  flow  model  are  grouped  into  a  single  matrix 
equation  of  the  type  (KARNI  et  al.  ,  1985) 

dP 

*  £P  =  Q  +  §.  (3) 

Here,  P  is  the  n-column  pressure  vector,  Z  is  the  symmetric  nxn  fluidity 
(inverse  of  resistivity)  matrix,  £  is  the  symmetric  nx«  compliance  matrix,  Q 
is  the  w-column  flux  vector,  and  S  -  the  n-column  source/sink  vector  which 
either  adds  or  drains  out  fluid  from  the  compartmental  domains.  If  5  =  0,  the 
system  is  conserved and  equation  (3)  reduces  to  the  conservation  of  mass, 
also  known  as  the  continuity  equation.  The  numerical  solution  presented  here 
is  for  a  conservative  system.  Solutions  for  non-conservati ve  systems  are  dis¬ 
cussed  elsewhere  (SOREK  et  al.,  1985). 

In  general,  equation  (3)  is  an  inhomogeneous,  ordinary  differential 
matrix  equation  of  the  first  rank  with  respect  to  time  (t).  It  is  also  non¬ 
linear,  since,  if  homeostasis  is  taken  into  account,  the  sensory  and  endo¬ 
crinological  biofeedback  control  mechanisms  turn  the  coefficient  matrices  c 
and  z  functions  of  the  dependent  vector  p  and  of  the  independent  scalar  t. 
Thus,  in  the  general  case, 
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£  =  g[_p(t),tj ;  z  =  z  \pjt),q 

There  are,  however,  some  passive  cases  that  can  be  approximated  by  the 
linear  problem,  i.e.  when  C  and  z  assume  constant  values.  Moreover,  even  if 
the  problem  is  non-linear  we  have  to  pass  through  the  linear  case  and  pro¬ 
ceed  to  the  non-linear  case  by  applying  incremental  perturbation  techniques. 
The  rest  of  the  discussion,  therefore,  focuses  on  the  linear  case. 

The  procedure  for  solving  the  linear  problem  is  as  follows: 

(a)  Given  information  about  fluxes  and  pressures,  equation  (3)  can  be  solved 
to  yield  the  resistances  and  compliances.  This  is  often  referred  to  as 
parameter  identification  or  model  calibration. 

(b)  When  the  elements  of  Z_  and  c_  have  been  evaluated,  the  pressure  waves 
P  =  P(t)  or  the  fluxes  Q  =  Q(t)  can  be  determined. 

(c)  With  the  pressures  and  the  compliances  solved,  information  about  the 
volume  changes  and,  under  certain  assumptions,  also  of  the  compartmental  dis¬ 
placements  can  be  obtained. 

We  shall  next  apply  this  to  the  case  of  a  seven-comoartmental  (n=7) 
model  that  represents  the  intracranial  cerebrovascular  fluid  system. 

2.  The  Seven-Compartmental  Model  with  Constant  Resistances  and  Compliances 

The  numerical  example  that  we  chose  to  calculate  is  based  on  AGARWAL'S 
six-compartmental  model  ( loc .  cit.,  1971),  for  which  no  numerical  results  were 
given,  with  the  addition  of  the  brain  tissue  (B)  compartment;  seven  altogether 

(Fig.  1). 


Fig.  1 
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The  lumped  resistances  are:  between  the  artery  and  capillary  compart¬ 
ments  the  capillary  and  cerebrospinal  fluid  compartments  (Rcp)>  the 

capillary  and  brain  tissue  compartments  (Rno),  the  capillary  and  vein  compart- 

Ld 

ments  ( RCy) ,  the  brain  tissue  and  vein  compartments  (RpvK  the  cerebrospinal 
fluid  and  brain  tissue  compartments  (RpBK  the  vein  and  venous  sinus  compart¬ 
ments  ^VSK  the  cerebrospinal  fluid  and  the  venous  sinus  compartments  (RpsK 
and  between  the  venous  sinus  and  the  jugular  bulb  compartments  (R.,  J ;  altogether 
nine  resistances.  In  the  figure,  the  capillary  compartment,  likewise  the  R.c 
resistance,  are  divided  into:  the  choroid  plexuses  -  those  tufts  of  small 
capillary  vessels  inside  each  of  the  four  ventricles,  and  the  capillary  system 
outside  the  ventricles.  However,  in  the  equations  to  follow,  only  the  combina¬ 
tions  in  parallel  of  the  resistances,  namely  (rac)  +  (rac^  irac)  and 

-1  -1  -1 

(R'cv)  +  -  <RcV*  appear,  so  that  the  combined  lumped  resistances 

Rac  and  Rcv,  into  and  out  of  the  capillaries,  suffice  for  our  purposes. 

The  resistances  R Rcp  and  RpB  are  identified  as  the  lumped  blood- 
brain  barrier;  the  lumped  blood-cerebrospinal  fluid  barrier,  and  the  lumped 
cerebrospinal  fluid-brain  barrier,  respectively.  Quantitative  studies  of 
these  barriers  indicate  that  on  the  gross  compartmental  level  they  assume 
large  values,  yet  they  cannot  be  regarded  infinite  (DAVSON,  1960). 

The  compliance  elements  C..  indicate  that  an  increase  in  volume  of  one 

id 

of  the  compartments  equals  the  volume  of  the  "cup"  formed  by  the  deformed  mem¬ 
brane.  This  volume,  in  turn,  equals  the  volume  displaced  from  the  neighbouring 
compartments,  all  within  the  rigid  container  of  the  skull  bones. 


First,  a  compliance  element  C  is  introduced  between  the  artery  and 
the  brain  tissue  compartments.  This  element  represents  the  overall  pulsatory 
effect  of  the  arteries  on  the  brain  tissue.  Next,  a  compliance  element  cCF 
is  inserted  between  the  choroid  plexuses  and  the  cerebrospinal  fluid  compart¬ 
ment.  The  CSF  system  and  the  brain  tissue  share  common  boundaries  -  at  the 
ventricles  and  along  the  subarachnoidal  space  -  which  are  deformable,  so  that 
a  compliance  element  CpB  is  introduced  between  them.  Finally,  as  result  of 
the  sharp  drop  in  pressure  along  the  cardiovascular  passage,  from  the  artery  to 
the  jugular  compartments,  additional  compliance  elements  CB]/  and  Cps  are  added 
between  the  brain  tissue  and  the  venous  compartments,  and  between  the  cerebro¬ 
spinal  fluid  and  the  venous  sinus  compartments,  respectively  (there  is  no  com¬ 
pliance  element  between  the  sinuses  and  the  large  jugular  veins  because  the 
pressure  there  is  very  small).  Altogether,  the  seven-compartmental  model  for 
the  cerebrospinal  fluid  system  described  above  assumes  five  compliance  ele¬ 
ments  of  the  type  .. 

The  mechanical  properties  of  'resistances’  and  'compliances'  are 
"symmetric"  with  respect  to  the  interchange  of  direction  between  one  compart¬ 
ment  and  its  neighbour.  This  is  the  outcome  of  the  law  of  action  and  reaction. 
In  formulae 


rac  ~  rca  rcf  rfc  eta' 


CAB  ~  °BA  CCF  =  CFC  etG’ 


The  chart  of  mean  pressure  variations  along  the  cerebrovascular  fluid 
conforms  with  the  pressure  profile  of  the  cardiovascular  system  cited  in  the 
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literature  (e.g.,  GUYTON,  1969,  Ch.  14).  In  the  arteries,  the  average  pres¬ 
sure  between  systole  and  diastole  is  100  mm  Hg.  It  drops  to  30  mm  Hg  in  the 
capillaries  including  the  choroid  plexuses;  10  mm  Hg  in  the  CSF  system,  9.5 
mm  Hg  in  the  brain  tissue,  9  mm  Hg  in  the  venous  system,  8  mm  Hg  in  the 
sinuses  and  1-2  mm  in  the  large  veins  -  jugular  and  spinal  -  leading  to  the 
vena  cava  (Fig.  1).  Volumes  of  the  compartments  are  also  recorded  in  the 
figure  as  much  as  they  are  documented  in  the  literature. 

For  the  flux  matrix  Q,  we  assume  the  following: 

The  heart  pumps  blood  into  the  artery  compartment  (arteries  and  arteri¬ 
oles)  at  the  rate  of  approximately  750  ml/min.  The  larger  amount  of  blood  flows 
into  the  capillaries  branching  outside  the  cerebral  ventricles,  while  a  small 
amount  reaches  the  choroid  plexuses  which  are  the  capillary  zones  inside  the 
ventricles.  No  definite  information  is  available  about  the  partition  ratio  be¬ 
tween  the  flows  through  the  capillary  section  and  through  the  choroid  plexuses. 
In  our  scheme,  the  ratio  of  A=250:l  has  been  postulated,  namely  747  ml/min  of 
blood  being  carried  into  the  capillaries  against  3  ml/min  entering  the  choroid 
plexuses  in  all  the  four  ventricles.  However,  as  mentioned  before,  only  the 
compound  resistances  RAC  and  Rcv\  so  also  the  fluxes  Qc,  Qy,  enter  the  calcula¬ 
tions  and  no  use  is  made  of  the  ratio  X  later  on. 

The  ultrafiltration  of  the  blood  at  the  choroid  plexuses  diverts  a  flow 
of  0.3  ml/min  to  the  CSF  compartment.  This  is  a  figure  extensively  quoted  in 
the  literature.  From  the  capillaries,  blood  is  recollected  into  the  venous 
system  (venules  and  veins).  A  minute  fraction  escapes  as  interstitial  fluid 
to  the  extracellular  region  of  the  brain  tissue  but  because  of  the  blood-brain 
barrier,  it  is  hardly  measurable  in  a  tenth  of  ml/min.  It  is  therefore  marked 
as  0.0  ml/min  in  Figure  I,  likewise  with  the  flow  from  the  cerebrospinal  fluid 


Li 

compartment  to  the  brain  tissue  compartment.  The  venous  compartment  also 
regains  the  3.0-0. 3=2. 7  ml/min  of  blood  which  does  not  escape  from  the 
choroid  plexuses  into  the  CSF  system,  thus  totalling,  with  the  flow  from  the 
CSF  system  into  the  brain  tissue,  749.7  ml/min  that  proceed  to  the  venous 
sinus  compartment.  The  drainage  of  the  CSF  into  the  sinuses  adds  the  0.3  ml/ 
min  which  closes  the  loop  flow  of  the  CSF  from  the  choroid  plexuses  to  the 
sinuses.  The  rejoining  of  the  CSF  drives  the  original  flow  of  750.0  ml/min 
back  to  the  heart  through  the  jugular  and  the  spinal  veins.  This  scheme  im¬ 
plicitly  assumes  no  gains  or  losses  in  the  flow  from  one  compartment  to  the 
other,  namely  the  source/sink  matrix  S  =  0. 


3.  Parameter  Identification  of  the  Seven-Compartmental  Model 

Given  the  above  data,  the  explicit  expressions  for  the  fluidity  z 
matrix  and  the  compliance  C  matrix  are: 


t  [zv]  Drr]  = 


AC 


-Z 


AC 


0  \?CP+7FPi’7'BV. 1  ~ZCB  ~Z 


FB 


' ZBV  0 


AC  '  CB  L  AC  CF 


[z_+z_+z 


cb+zcf\ 


-z 


CF 


-Z 


cv 


-z 


FB 


-Z, 


CF  \^CF+ZFS+ZFb\ 


-Z, 


-Z 


BV 


~ZCV  0  \^CV+ZVS+ZBv\  ~ZVS 


-Z 


FS 


JVS 


0 


0 


0 


\^VSfZFs\ 


I 
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g-M-fri]  - 


CAB  ~CAB 


~CAB  \?AB+CFB+CBv\  0 


JCF 


~CFB  ~CBV 


■CCF  0 


- FB 


'BV 


~C"-  ]?cf+cfs+cfb\  0 


CF 


' BV 


-°FS  0 


(5) 


-c 


FS 


JFS 


The  pressure  P  vector  has  the  following  elements 


P  =  {pa>pb>pc>pf>pv>ps] 


and  the  flux  Q  vector  reads: 


Q  =  {Qa,  0,  0,  0,  0,  -$.} 


The  parameter  identification  is  performed  in  two  stages: 

I.  The  steady  state.  For  constant  pressures,  the  matrix  equation  to  be  solved 
is  the  particular  integral  of  the  differential  matrix  equation  (3),  namely 


Z  P*  =  Q* 


(6) 


(7) 


(8) 
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Hence,  P*  and  <5*  indicate  the  average,  time- independent  compartmental  pres¬ 
sure  and  flux  vectors,  respectively,  the  values  of  which  are  given  in  Figure 

I.  For  the  linear  problem,  equation  (8)  is  a  solution  of  the  steady  state; 
it  is,  by  the  uniqueness  theorem,  also  the  solution. 

II.  The  non-steady  state.  Once  the  matrix  Z_  is  determined  from  the  steady 
state,  we  insert  it  into  equation  (3)  rearranged  to  read 

dP(t) 

g—ft —  =  Q(t)  -  Z  Pit)  f£,  Z  =  const.)  (3 ' ) 


The  solution  for  £  depends  on  the  availability  of  information  about  the  time 
dependency  of  both  the  flux  vector  Q(t)  and  of  the  pressure  vector  Pit).  In 
addition,  the  boundary  conditions  have  to  be  in  the  form  of  pressure  variation 
with  time  in  order  to  solve  uniquely  for  the  compliances. 

To  proceed  with  the  numerical  solution  of  these  equations,  we  replace 
the  partial  derivatives  by  implicit,  backwards  time  difference  quotients: 


dPAlt> 

dt 


5 1 


PA( ty+kht) ) 


(9) 


At  the  same  time,  we  express  all  the  space  derivatives  at  the  new  time, 
namely  at  n+1.  For  example 


PA(t)-PC(t) 


AC 


PnA+1(t)~flc1(t) 


R 


AC 


etc. 


The  differential  matrix  equation  now  turns  into  a  system  of  six  sets  of 
algebraic  equations,  each  having  the  form 
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a.cP^+b  .P+'+cP^+d.^+e  .Jl+1+f.p*'hl  = 

i  A  i  B  i  C  z  F  z  S  J  zV 


(i  =  1...6) 


(10) 


For  example,  for  the  artery  chamber  (i=l)3  we  have 


In  this  system  of  equations  we  have  pressure  values  at  times  n  and  n+l  (namely, 
at  times  tQ+nM  and  tQ  +  (n+l)kt  respectively),  values  of  resistances  H^.  which 
we  already  know  from  the  inverse  solution  of  the  steady  state,  and  the  values 
of  the  five  compliances  as  unknowns.  Given  the  values  of  pressures  at  times  u 
and  (n+l),  we  can  explicitly  solve  these  equations  for  the  five  compliances  cAli, 

CCFJ  CFB3  CFV*  and  CBV' 

Using  the  data  base  of  Figure  1,  the  compartmental  resistances  turn  out  to 
be: 


II 

Cj 

QC* 

0.0933 

mmHg/ml/min 

RCF  = 

66.67 

rr 

(lumped  blood-cerebrospinal  fluid 

barrier) 

hcb  = 

13,300.0 

tt 

(lumped  blood-brain  barrier) 

II 

0.0280 

tt 

RFB 

13.30 

tt 

(lumped  cerebrospinal  fluid-brain 

barrier) 

rfs 

7.62 

tt 

rbv 

12.  77 

tt 

Rvs 

0.0013 

tt 

RSJ 

0.0080 

tt 
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To  solve  numerically  for  the  compliances,  the  time-dependent  compart- 
mental  pressures  have  to  be  taken  from  clinical  measurements.  For  our  solu¬ 
tion,  we  chose  those  of  HAMIT  et  at.  (1965)  discussed  below.  Solving  equations 
(10)  leads  to  the  following  values: 

CAB  =  0.0012  ml/mmHg 
CCF  =  0.0357  » 

CBy  =  0.3750  " 

Cps  =  0.0494 

CpB  =  0.2090  " 

This  concludes  the  calibration  of  the  seven-compartmental  model  of  the 
cerebrovascular  fluid  system  with  constant  resistances  and  compliances. 


4.  The  Quasi -Steady  State  Pressure  Waves 

Once  the  linear  model  has  been  calibrated,  it  is  possible  to  obtain 
solutions  for  the  pressure  waves  in  all  the  compartments. 

The  solution  of  the  inhomogeneous  differential  matrix  equation  for  P 
in  incremental  form  is 


p(t+te)  =  erp 


-  +  Z -2q, 


?(-*■>  [l(t)  -  £  'S.(t)  |  '  =  *-{t) 


(ll) 


Here,  At  =  denotes  a  time  step  between  the  frontier  time  level  t^+j 


and  backtime  level  tfc.  The  matrix  k  equals 


k  =  At  •  C  Z 


(12) 


% 
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The  rational  expansion  of  exp  (-<_)  leads  to  the  following  formula 


I-(l-B) K 

exP(~ %>  3  TT'OiT^' 


(0z  0  <  1) 


in  which  I  is  the  nxn  unit  matrix.  A  substitution  of  equations  (12),  (13)  into 
equation  (11)  yields 


Z<Mt>  ‘  <L  *  <*''J  [i  [j 


=  -U). 


~  ±(t) 


The  coefficient  0  controls  the  type  of  solution  which  evolves  in  time. 

When  0=0,  we  have  an  explicit  scheme;  0=1,  an  implicit  scheme,  and  0  <  0  ■  1  is 
the  mixed  scheme.  Thus,  with  the  choice  of  0,  the  pressure  waves  in  the  various 
compartments  can  be  calculated  from  equation  (14). 

Yet  there  is  another  state  between  the  steady  state  and  the  non-steady 

state.  If  we  consider  the  contribution  of  the  compliance  term  to  be  negligible, 
dP 

C  ~  «  0  (which,  in  effect,  is  the  case  once  we  insert  the  numerical  values 
*  at 

of  Cjj  into  equation  (31)),  we  obtain  the  relation 

Z  P(t)  =  Q(t)  (15) 


Equation  (15)  "looks"  like  equation  (8)  -  the  condition  for  the  steady  state 
flow  -  but,  this  time,  neither  the  pressure  P  nor  the  flux  U  assume  constant 
values;  rather,  each  of  them  is  a  function  of  time.  This  is  the  guae.i-e.teaJy 
state  and  the  pressure  waves  can  be  calculated  from  a  simplified  formula, 

namely 
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ttuet)  *  k\t)  (14'> 

HAMIT  et  al.  (loo.  cit.)  performed  simultaneous  recordings  of  EC6,  PCG, 
arterial ,  brain,  cisternal  and  venous  (strictly,  venous-sinus)  pressure  waves 
in  anesthesized  dogs  (Fig.  2b).  The  recorded  waves  were  in  the  6-8  cycle-per- 

Fig.  2 

minute  (c.p.m. )  range  upon  which  the  faster  cardiac  waves  -  also  termed  "pulse 
waves"  -  were  superimposed.  In  the  classification  of  the  intracranial  pressure 
waves,  the  slow  ones  correspond  to  the  Lundberg  B-waves  (LUNDBERG  ct  al.,  loo. 
oit.,  1974)  and  the  faster  pulse  waves  to  the  Lundberg  A-waves. 

The  arterial  B-wave  shows  an  almost  ideal  sine  wave  pattern  varying  be¬ 
tween  110-140  mmHg.  We  took  it  as  the  excitation  pressure  wave  in  the  artery 
compartment  at  frequency  0.144  Hz  or  a  7  second  period,  and  calculated, 

by  means  of  equation  (14'),  the  rest  of  the  compartmental  pressure  waves.  The 
computer  plotting  of  the  venous-sinus  pressure  wave  Pq(t)  shows  a  sine  pattern 
of  the  same  frequency  at  amplitudes  varying  between  2.8  -  7.2  mmHg  (Fig.  2a) 
which  is  in  excellent  agreement  with  the  Pn  -  wave  measured  by  HAMIT  at  al. 

The  compartmental  analysis  also  allows  estimates  about  some  kinematical 
changes  which  take  place  with  the  pressure  waves.  On  the  lumped-compartmental 
scale,  the  volumetric  deformation  of  the  brain  tissue  (/>')  compartment  is  re¬ 
presentable  by  the  compliance  coefficient  CFB-  As  listed  above, 

C  =  0.2090  ml /mini!  ij 

F  B 


zz 


The  pressure  difference  (c»/.Fig.  1)  ppg  =  Pp-iJB  =  0.5  mmHg.  Following 
HAKIM  et  al.  (loo.  ait .,  1976)  we  assume  the  brain  tissue  to  lie  in  a  spheri¬ 
cal  layer  between  the  inner  radius  of  the  ventricles  r.  and  the  outer  radius 

1  v 

of  the  subarachnoidal  space  r  ;  the  ratio  of  the  radii  to  be  =  4,  and  the 

i 

volume  Vq  =  600/2i r  mi.  Thus,  the  radial  displacement  of  the  ventricles  is 

found  to  be 

-4 

Ar •  =  0. 128  x  10  mm 

HEIFETZ  and  WEISS  (1981)  have  shown  in  two  patients  that  after  raising 
the  cerebrospinal  fluid  pressure  by  15-20  mmHg  (namely,  increasing  30-40 
times  the  value  of  PrD  cited  before),  measurable  changes  occurred  on  electri- 
cal  capacitance  strain  gauges  fastened  to  the  skull.  Converting  their  results 
to  displacements,  the  movements  of  the  skull  in  response  to  the  pressure  ele¬ 
vation  was  0.00078  mm  in  one  case  and  0.00372  mm  in  the  other.  Multiplying 
the  above  calculated  value  of  Ar^  and  ppg  =  0.5  mmHg  by  the  factor  30-40,  we 
find  a  very  good  agreement  with  the  measurements  of  HEIFETZ  and  WEISS. 

A  lateral  skull  motion  of  the  order  of  a  few  microns  was  also  recorded 
by  IVAN  et  al.  (1983)  using  electrical  resistance,  high  extension  rubber 
strain  gauges  (Peekel  type  20S)  placed  over  the  skull  sutures.  Earlier, 

FRYMAN  (1971),  who  used  spring  dial  gauges,  recorded  lateral  motion  of  the 
temporal  bones  in  the  order  of  10-18  microns  (p).  ALLEN  at  al.  (1983),  who 
analyzed  cine-CT  scans,  reported  slow  rhythmical  deformations  of  the  ventri¬ 
cles  -  the  third  ventricle  in  particular  -  of  the  order  of  0.1  -  1.0  p,  in 
conformity  with  the  above  data  at  least  to  the  order  of  magnitude. 


The  compartmental  approach  is  a  useful  tool  in  the  modelling  of  intra¬ 
cranial  fluid  dynamics  as  far  as  the  time  dependency  is  concerned.  Its  major 


23 


/ 


drawback  is  that  it  does  not  relate  events  to  their  spatial  configuration 
since,  by  definition,  the  lumping  of  the  parameters  is  space-independent. 

For  the  space-time  modelling  of  brain  tissue  mechanics,  we  have  to  revert  to 
the  continuum  or  distributed  parameter  modelling  which,  on  the  macroscale, 
considers  the  brain  tissue  single  phased  possessing  the  averaged  property  of 
viscoelasticity,  and  on  the  microscale,  analyzes  it  as  a  multiphasic  system 
-  three  at  least  -  of  neurons,  glia  and  interstitial  fluid.  Here,  much  more 
data  is  needed  to  calibrate  the  models  and  before  this  is  reached,  the  con¬ 
tinuum  modelling  will  have  to  stall. 

Finally,  the  illustrated  example  of  the  compartmental  modelling  here 
is  seven-compartmental .  If,  in  view  of  additional  data,  a  model  of  more  com¬ 
partments  is  favoured;  or  the  choice  of  resistances  and  compliances  altered, 
nothing  changes  in  the  methodology  described  above  or  in  the  computer  pro¬ 
gramme  which,  anyway,  was  programmed  for  nxn  matrices.  The  authors  will  be 
grateful  if  information  about  other  data  bases  be  brought  to  their  attention. 
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Fig.  2.  a.  Computed  Ps  wave  from  a  PA  excitation  wave 
(0.144  Hst  110-1-10  mm  Hg). 
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Abstract  -  A  lumped  parameter  compartmental  model  for  the 
cerebrovascular  fluid  system  is  constructed  and  solved  for  the 
general  linear  problem  of  a  nonsteady  flow  with  constant  re¬ 
sistances  and  compliances.  The  model  predicts  the  intracranial 
pressure  waves  in  the  various  compartments  of  the  brain  in 
response  to  pressure  changes  in  the  vascular  system. 

NOMENCLATURE 


(*) 

A  '  sensitivity  matrix 

A  arterial 

B  deformationl  flow 

B  brain 

C  compliance  matrix 

C  capillary 

0  difference  vector 

F  cerebrospinal  fluid  (CSF) 

J  jugular  bulb 

l’  pressure 

J  flux  vector 

H  resistance 

<j  venous  sinus 

*In  print,  letters  underlined  once 
♦♦Deceased 


•n 

(  )  transpose 

L  time 

v  venous,  volume 

Z  fluidity  matrix 

u  ratio  of  CSF-brain  resistance 

to  vein-venous  sinus  resistance 

3  ratio  of  CSF-brain  resistance 

to  capillary-brain  resistance 

y  ratio  of  capillary-vein  res¬ 
istance  to  brain-vein  resistance 

0  coefficient 

II  n-set  pressure  matrix 

or  twice  -  boldface  type 
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INTRODUCTION 

This  paper  summarizes  the  first  stage  of  a  research  aimed  at  modelling  the 
movement  of  fluids  and  chemicals  in  the  cerebrovascular  system.  Here,  we 
present  the  solution  of  the  non-steady  flow  with  constant  resistances  and 
compliances;  namely,  the  general  linear  problem  -  in  its  entirety. 

The  lumped- parameter  compartmental  model  of  the  cerebrovascular  system 
is  the  first  step  towards  a  more  comprehensive  modelling  of  the  cerebral 
content  considered  as  a  material  continuum  both  on  the  macro  and  on  the 
micro-scales.  By  the  compartmental  approach,  the  intracranial  content  is 
divided  into  a  number  of  subunit  compartments  the  behavior  of  each  of  which 
is  represented  by  a  single  pressure  parameter  and  by  a  single  fluid  dis¬ 
charge,  both  functions  of  time  but  not  of  space.  The  resistance  to  flow 
due  to  a  particular  vessel  type  is  lumped  at  the  outflow  of  its  compartment. 
Likewise,  the  intergrated  change  in  volume  of  each  compartment  of  the 
system  is  representable  as  an  overall  compartmental  property  and  the 
functional  interaction  between  the  components  of  the  lumped  parameter 
system  is  assumed  to  be  at  the  interfaces  between  adjacent  compartments. 

Monro's  (1783)  first  model  of  the  intracranial  cavity  was  bi -compartmental 
brain  and  blood  as  two  almost-incompressible  phases.  Kellie  (1824), 

Monro's  pupil,  modified  this  hypothesis  by  assuming  three  instead  of  two 
material  compartments,  namely  arteries,  veins  and  brain  tissue.  The  Monro- 
Kellie  doctrine  prevailed  to  this  century  and  was  only  relaxed  in  stages. 

The  number  of  fluid  compartments  was  increased  to  six:  arterial,  capillary, 
venous,  venous  sinus,  jugular  bulb  and  cerebrospinal  (CSF)  (Agarwal,  1971). 
Yet,  the  fluid  itself  remained  incompressible.  More  recent  approaches  relaxed 
the  latter  assumption  and  the  fluid  was  taken  as  linearly  compressible, 
namely  the  change  in  pressure  and  the  change  in  volume  are  proportional  to 


each  other,  and  the  coefficient  of  proportionality,  the  bulk  modulus  or  its 
inverse  -  the  compliance,  are  constants.  As  to  the  brain  tissue,  in  essence 
a  multiphasic  material  continuum,  the  experimental  results  pointed  at  non¬ 
linear  elasticity  of  the  tissue  or,  alternatively,  the  nonlinearity  of  the 
compliance  (Marmarou  et  aZ.,1975;  Miller, 1975;  Hakim  et  at., 1976;  Lewer 
Allen  and  Bunt,  1978;  Chopp  and  Portnoy,  1980).To  overcome  the  "non-linearity" 
of  a  single  coefficient,  more  complex  models  of  single-phased,  multi¬ 
parameter  viscoelastic  materials  were  introduced  such  as  the  one  consisting 
of  four  viscoelastic  coefficients  known  as  the  "three-parameter  solid", 
coupled  dynamically  with  another  elastic  element  (Pamidi  and  Advani ,  1978). 

Our  model  consists  of  seven  compartments , namely  the  six  compartments 
listed  above  and  the  brain  tissue  compartment  (Fig.l).  The  lumped 
resistances  ?.re:between  the  Artery  [a)  and 

Fig.  1 

Capillary  (c)  compartments  (f^),  the  Capillary  ( C )  and  Cerebrospinal 
Fluid  (f)  compartments  (f^),  the  Capillary  (c’)  and  Brain  Tissue  {/<) 
compartments  (RCB) ,  the  Capillary  (c)  and  Vein  (v)  compartments  (A’(?y),  the 
Brain  tissue  («)  and  Vein  (v)  compartments  (i<BV).  The  Cerebrospinal  Fluid 
(F)  and  the  Brain  (B)  compartments  ,  the  Vein  ('/)  and  Venous  Sinus  (:<') 

compartments  (/?.„),  the  Cerebrospinal  Fluid  (f)  and  the  Venous  Sinus  (F) 
compartments  (FfS),  and  between  the  Venous  Sinus  (f)  and  Jugular  Bulb  (<J) 
compartments  (z?SJ);  altogether  nine  resistances.  In  the  Figure,  the 
Capillary  compartment,  likewise  the  f  resistance,  are  subdivided  into 
the  chambers  of:  the  choroid  plexuses  -  those  tufts  of  small  capillary 
vessels  inside  each  of  the  four  cerebral  ventricles,  *nd  the  capillary  system 
outside  the  ventricles.  However,  in  the  equations  to  follow,  only  the 
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combinations  of  flj"3  +  /rj~2=  r'J,  and  r'^  +  ff"“ J  =  ff‘3  appear  so  that  the 
lumped  resistances  R^c  and  r  ,  into  and  out  of  the  capillary  com¬ 
partment,  prevail.  Nevertheless,  it  is  justified  to  conserve  the  individual 
components  of  R^  and  R because  the  latter  relate  to  anatomically 
distinct  features  through  which  different  flows  take  place.  We  note  that 
the  lumping  is  a  consequence  of  the  fact  that  the  pressure  difference 
across  the  components  is  the  same. 

The  resistances  Rcg  ,  Rcp  and  Rpg  are  identified  as  the  blood-brain 
barrier;  the  blood-cerebrospinal  fluid  barrier  and  the  cerebrospinal 
fluid-brain  barrier,  respectively. 

The  compliance  elements  C..  indicate  that  an  increase  in  volume  of  one 

i,j 

compartment  equals  the  volume  of  the  cup  formed  by  the  deformed  membrane. 
This  volume,  in  turn,  equals  the  volume  displaced  from  the  neighboring 
compartments,  all  this  within  the  overall  rigid  container  of  the  skull 
(the  Monro-Kellie  doctrine). 

In  the  non-steady  state,  the  deformabil ity  of  the  compartments  due  to 
the  pulsatory  motion  of  the  arteries  is  taken  into  account.  Thus,  first, 
we  introduce  a  compliance  element  C  between  the  artery  and  the  brain 
tissue  compartments.  Next,  the  capillary  system  is  considered  non-deformable 
so  we  do  not  insert  compliances  between  this  compartment  and  any  of  its 
neighbors.  The  choroid  plexuses,  however,  although  being  capillary  in 
character,  possess  other  material  properties  and  can  -  and  in  fact,  do  - 
convey  pulsations  to  the  CSF  system  (Bering  (1955))  so  that  a  compliance 

element  is  introduced  between  them.  Further,  the  CSF  system  and  the 

■  7- 

brain  tissue  share  common  boundaries,  at  the  ventricles  and  along  the  sub¬ 
arachnoidal  space,  which  are  mechanically  deformable.  A  compliance  element 
is,  therefore,  inserted  between  the  two.  Finally,  additional 
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compliance  elements  CgV  and  t'FS  are  introduced,  respecti vely ,  between  the 

brain  tissue  and  venous  compartments,  and  between  the  CSF  and  venous  sinus 
compartments  to  account  for  possible  changes  in  fluid  accumulation 

between  the  considered  compartments,  and  in  order  to  reduce  the  high 
pressure  in  the  cardiovascular  passage.  There  is  no  compliance  element 
between  the  sinuses  and  the  large  jugular  veins  because  the  jugular 
bulb  is  assumed  rigid,  with  very  small  storage  changes.  Altogether, 
in  our  presentation,  we  assume  five  compliance  elements  between  the 
cerebrovascular  elements. 

The  mechanical  properties  of  ‘resistances'  and  'compliances'  are 
"symmetric"  with  respect  to  the  interchange  of  direction  between  one 
compartment  and  another.  These  are  outcomes  of  the  law  of  action  and 
reaction.  In  general  formulae, 
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Three  results,  obtained  by  running  the  model  presented  here  for  cases 
of  interest,  indicate  that  the  model,  even  in  its  present  form,  can 
simulate  certain  features  actually  observed  clinically.  The  first 
relates  to  the  observed  displacement  of  the  main  ventricle  peripheral 
boundary.  By  computing  the  value  of  the  compliance  between  the  F-compartment 
and  the  Brain  Tissue  one,  approximating  the  configuration  of  the  F- 
compartment  to  a  sphere,  Kami  et  al.  (1986)  showed  that  this  displacement 
is  of  the  order  of  magnitude  actually  observed  by  Lewer  and  Bunt  (1978). 
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The  second  indication  relates  to  the  pressure  variations  in  the 
jugular  bulb.  Sorek  et  al.  (1986a),  by  employing  the  model  presented 
here,  predict  that  a  negative  pressure  appears  during  certain  portions 
of  the  pressure  cycle  in  the  jugular  bulb  compartment.  Such  negative 
pressure  is  indeed  observed  in  reality.  Also,  the  model  predicts, 
during  certain  portions  of  the  cycle,  zero  pressure  in  the  sinus,  a  fact 
that  is  also  accounted  for  the  action  of  the  heart. 

Finally,  as  shown  in  this  paper,  the  model  predicts  the  excessive 
accumulation  of  CSF  in  the  F-compartment  (known  as  Normal  Pressure 
Hydrocephalus),  as  a  result  of  clogging  of  the  passage  between  the 
capillaries  and  the  Brain  Tissue  compartments .  Sorek  et  al.  (1986c) 
report  observing  this  relation  and  discuss  the  use  of  it  in  relieving 
the  pressure  by  transplanting  a  shunt. 


COMPARTMENTAL  MODEL  EQUATIONS 


The  assumptions  underlying  the  hydrodynamical  model  of  the  circulatory 
cerebrovascular  system  may  now  be  summarized  as  follows: 

(1)  Functional  brain  elements,  represented  by  compartments ,  are  character¬ 
ized  by  'mean'  or  'lumped  parameter'  quantities,  of  interest  to  us  here 
are  pressure  and  fluid  flux. 

(ii)  The  fluid  is  single  phased,  incompressible  and  Newtonian.  From  the 
physical  point  of  view,  there  is  no  differentiation  between  blood  and  CSF. 


(iii)  The  flow  is  laminar  and  the  relation  between  the  flow  rate  and 
pressure  drop  is  linear,  i.e. 


o  = 


A  r 

if 


(2) 
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where  R  is  the  resistance. 

(iv)  For  the  distensibi 1 ity  (elasticity)  of  the  element  of  brain,  a  linear 
function  holds  between  the  volume  change  AP  and  the  pressure  drop  AP 
namely 

AV  =  CmAI'  (3) 

where  the  coefficient  of  proportionality,  c.  ,  is  the  compliance. 

(v)  The  system  is  isothermal  and  thermodynamically  stable. 

The  explicit  continuity  equations  for  the  balance  of  mass  for  the 
various  compartments  now  read: 

For  the  Artery  ( a )  compartment; 


(7) 


For  the  Vein  (/)  compartment: 
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P  -P  P  -P  p  -P 

S±L  +  =  tvj’s  d__  (y  _p 

Rcv  RBV  Hvs  Vli  dt  v  R 


(8) 


For  the  Venous  Sinus  (5)  compartment: 


Fv~Fs  RF~RB  PB~PT  1 

- —  ^  _L tl.  =  - £  +  p  d  (p  _p  ) 

RVS  RFS  Rp,r  RR  dt  P  F 


(9) 


For  the  Jugular  Bulb  (J)  compartment: 


FS~P'J 
~RSJ  '' 


(10) 


Adding  (4  )  through  (10),  making  use  of  (  1  ),  yields 


qa-qj  (11) 

which  is  a  remanifestation  of  the  Monro-Kellie  doctrine. 

Equations  (4)-(10)  can  be  grouped  into  a  single  matrix  equation  of  the 
form 

dP 

£  ft  +  l  £  =  q  (12) 

where  P  is  the  pressure  column  vector: 


1‘  {PA 


P  P  P 

B>  C*  Fy 


(13) 
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Z  is  the  fluidity  (inverse  of  resistivity)  matrix 


ts: 

1 

0 

~ZAC 

0 

0 

0 

0 

t ZCB+ZFB+ZBV ^ 

~ZCB 

~ZFR 

-ZRV 
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-he 
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- ZBV 
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0 
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- ZFS 

~ZVB 

[z^T+zf,J 

(14) 


C  -  the  symmetric  compliance  matrix 
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(15) 


and  Q  is  the  flux  column  vector 


Q  =  {Qa  >  0  ,  0  ,  0  ,  0  ,  -Q.)  (16) 

The  resistance  does  not  appear  in  equation  (14)  as  it  is  readily 
obtainable  from  equation  (10). 
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The  matrix  equation  (12)  comprises  nine  resistances  and  five  compliances 
listed  in  equation  (1).  Even  for  the  steady  state,  whereby  all  the 
derivatives  with  respect  to  time  vanish  identically,  and  equal  ion  (12) 
reduces  to 

Z  P  =  Q  (17) 

we  are  left  with  the  eight  resistances  (as  Rc  is  determined  separately) 

bj 

against  six  independent  balance  conditions  for  the  various  compartments  of 
the  cerebrovascular  model,  excluding  the  jugular  bulb.  However,  by  virtue 
of  the  Monro-Kellie  doctrine  which  assumes  an  absolute  rigidity  of  the 
cranial  vault,  actually  only  five  equations  of  the  six  are  independent. 

Thus,  the  redundancy  of  the  system  (17)  is  three,  and  three  additional 
conditions  are  needed  to  solve  the  set.  One  of  them  is  the  flux  Qp-  0.3  ml/min 
which  can  be  taken  as  pivot  value  with  high  credibility  from  the  literature 
(Cutalar,  1968).  Two  more  conditions  have  to  be  stipulated  by  physiological 
data.  They,  however,  do  not  assume  specific  values;  instead  they  are  known 
to  be  within  certain  limits.  We  shall  have,  therefore,  to  accompany  their 
determination  by  a  sensitivity  analysis,  sweeping  the  entire  range  of 
variation  of  possible  values,  as  discussed  in  the  following  Section. 

EVALUATION  OF  THE  FLUIDITY  MATRIX 

The  intracranial,  compartmental  flow  problem  of  the  cerebrovascular 
system  lies  in  the  solution  of  the  set  (12)  in  its  entirety.  This  is  an 
inhomogeneous,  first-order,  ordinary  differential  matrix  equation  with 
respect  to  time.  In  the  previous  Section,  we  chose  to  refer  to  the 
matrices  £  and  Z  as  constant,  a  tact  that  made  (12)  a  set  of  linear 
equations.  The  respective  mean  values  were  taken  as  these  constants, 
as  a  first  approximation.  However,  in  principle,  the  coefficient  matrices 
£  and  Z  could  be  functions  of  the  dependent  vector  parameter  P  which, 


in  turn  depend  on  the  time,  t  .  ihus,  basically  the  problem  is  non-l inear. 
There  are,  however,  some  passive  cases  which  can  be  approximated  by  the 
linear  problem  namely,  when  £  and  Z  assume  constant  values.  At  any  event, 
even  while  dealing  with  the  non-linear  problem,  we  shall  have  to  pass 
through  the  linear  case.  To  the  latter,  therefore,  our  attention  is  now 
focused. 

The  procedure  of  solution  of  the  linear  case  is  as  follows: 

(a)  Given  information  about  fluxes  and  pressures,  equation  (12)  can  be 
solved  to  yield  values  for  resistances  and  compliances.  This  is  often 
referred  to  as  parameter  identification  or  model  calibration. 

(b)  When  the  R.. s  and  C.. s  have  been  evaluated,  the  changes  in 

t'tf  I'd 

pressures  and  rates  of  fluid  flow  can  be  determined  for  variations  in  the 
rate  of  fluid  discharge  from  the  heart  pump  and  from  the  venous  outlet. 

(c)  With  the  pressures  and  compliances  already  solved,  information  about 
the  volume  changes  and,  under  certain  assumption,  also  of  the  displacements 
of  the  compartments  can  be  obtained. 

Here,  we  confine  the  attention  to  the  parameter  identification  and 
proceed  in  two  stages: 

!•  S’fce  steady  state.  For  constant  pressures,  the  governing  equation  is 
the  particular  integral  of  (12).  The  matrix  equation  to  be  solved  is  (17). 

The  chart  of  mean  pressure  variations  along  the  cerebrovascular  fluid 
system  conforms  with  the  pressure  profile  of  the  cardiovascular  system 
cited  in  the  literature  (e.g.,  in  Guyton, 1969,  Ch.14).  In  the  arteries, 
the  average  pressure  -  between  systole  and  diastole  -  is  100  nin  Hg.  It 
drops  to  20  mm  Hg  in  the  capillaries  including  the  choroid  plexuses; 

10  mm  Hg  in  the  CSF  system,  9.5  mm  Hg  in  the  brain  tissue,  5  mm  Hg  in  the 
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venous  system, 2  mm  Hg  in  the  sinuses  and  1-2  mm  Hg  in  the  larger  veins  - 
jugular  and  spinal  -  leading  into  the  vena  cava  (Fig.l).  Volumes  of  the 
compartments  are  also  recorded  in  the  Figure  as  much  as  they  are  documented 
in  the  literature. 

Given  the  compartmental  mean  pressures  P*  (i=A,B,(',F,V,S,J ) ,  also 

introducing  the  abbreviations  P*  ■  for  the  (liltercmt'  . . .  nr> 

tj 

between  two  adjacent  interacting  compartments,  l'\j  »'  larly, 

~  ®i~^j  rt>presents  the  difference  between  the  mean  fluxes  in  com- 

A  A  A 

partments  i  and  j.  Because  Q„,  Q  and  Q,  are  a-priori  known,  and 

ft  j  ' 

the  mean  pressure  values  P^,  P^  and  P^ f  (as  all  other  mean  pressures) 
are  also  known,  the  following  elements  of  t  are  immediately  solved: 


(18) 


This  leaves  us  with  six  elements  of  A  against  five  equations  for  the 
remaining  compartments,  four  of  which  are  independent. 

To  account  for  the  double  redundancy,  we  introduce  the  scalar  coefficients 


“vs 


in  0) 


(19) 


RFB  _  7,_cb 
“ Cl)  ZFB 


i  <>  _  !'•  *>) 


(20) 


Here,  a  indicates  the  ratio  of  the  cerebrospinal  fluid-brain  barrier  to  the 
vein-venous  sinus  resistance,  whereas  a  is  the  ratio  of  the  cerebrospinal 
fluid-brain  barrier  to  the  blood-brain  barrier. 
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Let  us  now  examine  the  following  cases: 


(i)  0  <  8  <  »  . 

leads  to: 


Based  on  equation  (20),  solution  of  equation  (17) 


(21.1) 


ZFB  a P*vs2-P*b  Q1f 


7  =  a  n* 

JVS  aP  -P*  hap 
KS  FB 


ZFS  ^VS^FH  *  *FS  (ClPvS  Q*F~P*FB  Q A/[) 


Z  -  <s-pfb-Kb  .  ,  . 

CV  ^vS-PpB  QaF 


cv 


(21.2) 

(21.3) 

(21.4) 

(21.5) 


ZCB 


B 

a P*  -P 
VS  FB 


AF 


ZBV  = 


m 


a P 


VS  FB 


FB 

?BV 


0 * 
*AF 


(21.6) 

(21.7) 


The  lower  extreme  value  of  8=0  is  noteworthy  of  mention.  From  the 
definition  of  8  (Eq.  (20)),  it  can  assume  a  zero  value  in  two  cases. 

First,  pFB=0,  indicating  a  complete  rupture  of  the  cerebrospinal -brain 
barrier.  Here,  a =0  also  prevails  and  the  only  consistent  condition  for 
equation (21.3)  to  exist  is  pF=pB*  which  means  that  there  is  hydrodynamical 
equilibrium  and  absolute  stagnancy  of  flow  between  the  CSF  compartment  and 
the  brain  tissue  compartment.  To  avoid  this,  the  pressures  at  the  CSF 
compartment  and  the  brain  tissue  compartment  ought  to  differ  by  some 
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amount,  P*  ft  P*  .  In  the  literature,  it  is  customary  to  assume  that 
both  are  in  the  order  of  10  mm  Hg.  In  our  calculaticns,  we  arbitrarily 
used  the  mean  values  of  P^=10  nun  Hg  and  pg=  9.5  mm  Hg  (Fig.l)  to  enable 
flow  from  the  F  to  the  B  compartment. 

The  second  possibility  for  the  vanishing  of  6  is  -*■  °“  . 
Physiologically,  this  means  that  there  is  a  complete  blockage  of  the  blood- 
brain  barrier;  there  is  no  transmission  of  blood  nutritions  from  the 
capillaries  to  the  interstitial  fluid  of  the  brain  tissue,  in  short  - 
a  total  collapse  of  the  brain.  We,  therefore,  rule  this  case  out. 

There  is,  however,  one  more  case  of  an  admissible  non-singular  solution 
of  the  hydrodynamical  problem  for  tt=0  .  It  occurs  when  2 '.„j=l/RpS=0 
or,  alternatively,  from  equation (21 .4), when  n= ^a'PFi/^f'PVS  ’  This 
(if  Pp  >Pfi)  indicates  that  there  is  no  direct  drainage  of  the  CSF  into  the 
sinuses  through  the  arachnoidal  granulations  (the  pacchionian  bodies),  the 
CSF  compartment  considerably  increases  its  size  and  the  symptoms  are  those 
of  hydrocephalus. 

There  is  little,  almost  no  data  about  the  physiological  or  pathological 
ranges  of  a  and  B.  By  trial  and  error,  however,  we  could  assume  B  to 
be  in  the  range  of  1/1000.  This  value  led  to  results  that  have  actually  been 
observed  in  the  clinic.  At  the  same  time,  the  values  of  the  compliances  and 
resistances  obtained  were  found  to  be  insensitive  to  the  value  of  B,  (see 
Sorek  et  al.  1986b).  To  convert  this  into  corresponding  ratios  of  fluxes, 
designated  by  a  bar,  instead  of  ratios  of  resistances,  we  have  by  definition. 


so  that 


B  i 


*CB 
? FB 


a 


-  _  PVSZVS  _  PV3_  RJB  _  P_VS_  a 
PFBZFB  PFB  RVS  PPB 

P  7  P  R  V 

p  _  CB  CB  ‘jOB_  _FB_  _  CB 

P  ~  P  Z  P  R  ~  P 

FB  L  FB  *  FB  CB  FB 


(22) 


(23) 
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For  our  choice  of  B  =  1/1000,  B  =  0.021  (c/.Fig.l).  As  illustrations. 
Figures  2  and  3  present  the  computer  solution  of  the  set  (21)  for 
RnJa)  and  R„Ja)  respectively  in  the  particular  case  of  B  =  0.021 

Ld  to 

Fig.  2,  3 


(B  =  1/1000).  The  functional  dependence  of  1<CB  on  u  is  linear, 

whereas  that  of  R„.(a)  is  hyperbolic. 
bS 

A  general  sensitivity  analysis  for  the  variation  of  the  resistances 
R..  in  the  present  case  (i)  of  0  <_  B  <  00  may  follow  the  same  procedure 

ij 

as  for  case  (ii)  (B  -*•  °°)  discussed  below.  This  was  not  carried  out 
numerically  because,  at  a  later  stage,  due  to  the  vagueness  of  the 
assumptions  about  B  ,  the  complete  solution  for  the  resistances  R^j(a.  3/, 
rather  than  the  limited  solution  of  for  the  shove  values  of  B, 

t'J 

has  been  obtained.  The  three-dimensional  plottings  for  the  compartmental 
resistances  R..(a,$)-,  later  also  those  for  the  compliances  C..(a,  P>), 
will  be  reported  separately. 

We  now  turn  to  case: 

(ii)  B  -*■  »  .  In  addition  to  equation  (18),  the  remaining  elements  of  the 
matrix  Z  read 


ZFB  °  ;  ZF'J  P*s  J' 


Z  -  — 

vs  ‘  ife 


ZCV  ~  p*  +  V p*  " AF 
CV  Y  BV 


—  0 * 
3* 


7jbv  yZcv 


7  -  1  Q* 

CB  ~  p *  +  vp*  p*  *AF 

CV  yBV  “CH 


(24) 
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Here,  for  the  sake  of  a  smoother  solution,  a  new  parameter  y  instead  of 
a  has  been  introduced  as  the  ratio  of  the  resistances  and  rd„  , 

t  V  D  V 

namely 


Y 


_  RCV 


R 


BV 


(25) 


Strictly,  the  exact  value  of  the  ratio  y  ,  as  much  as  those  of  a  and 
0  ,  is  not  known.  However,  we  can  still  pursue  with  the  sensitivity 
assessment  of  the  z  elements  against  a  variation  in  y,  >'*,  and  P*  . 
Let  A  denote  a  deviation  from  the  mean  value,  then 


A?*  =  A P* 

CV  “  BV 


+  A P 


CB 


(26) 


It  follows  from  equations  (24  ), (25)  that 

°-Z  ■  &°-F 

where 

-Z  =  ^CV  3  A lCB  3  Ll‘uv^ 
D-f  =  *  *CB  3  ^ 


(27) 


(28) 


are  the  difference  vectors  of  the  Z  elements  and  of  their  arguments 
respectively,  and  a  is  the  square  matrix 
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BV 
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-(1+  0 
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AF 


(PCV+yPBV)  ‘ 


BV 

P* 
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CV 


BV 

* 


CV 
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+  Y 


CB  CB 


BV 
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f  Vy 


(29) 
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The  extreme  case  of  6  -*■  »  corresponds  to  a  complete  rupture  of  the 
blood-brain  barrier  (RCB=0)  causing  haemorrhage  to  the  entire  brain 
tissue.  This  is  bound  to  change  the  compartmental  pressure  distribution 
contrary  to  our  assumption  of  a  steady  state.  Use  of  equations  (27), 

(28), (29)  is.  therefore,  confined  to  the  onset  of  the  occurrence  only. 

II.  Non  steady  flow.  An  interim  solution  of  the  matrix  equation  (12) 
for  P=P(t)  and  Q-Q(t)  when  the  effect  of  C_  is  negligible  but  not 
zero,  has  been  discussed  elsewhere  (Kami  et  al.3  1986).  As  an  example, 
for  a  given  arterial  excitation  wave  /’ =/’ (t)  ,  the  venous-sinus  wave 
Ps(t)  was  calculated  and  checked  against  the  measured  Ps(t)  wave 
recorded  by  Hamit  et  al.,  1965  ,  showing  a  good  matching  in  wave  form, 
amplitude  and  frequency.  In  this  way,  after  the  model  is  calibrated 
and  the  resistances  R..  solved,  the  complete  evaluation  of  the  pulse 
wave  vector  P(t)  can  be  obtained. 

EVALUATION  OF  THE  COMPLIANCE  MATRIX 

We  now  proceed  to  the  complete  solution  of  the  non-steady  flow  problem 
for  the  compartmental  model  with  constant  coefficients,  namely  the  general 
linear  problem. 

Equation  (12)  is  rearranged  to  read 

dP(  t) 

C  —  =  Q(t)  -  Z  P(t)  =  B(t)  (C,  Z  =  const.)  (30) 

Information  scarcely  exists  about  the  time  dependency  of  the  vector  q(t) 
at  all  compartments.  We,  therefore,  begin  with  the  average  pressure 
vector  P*  and  the  average  flux  vector  Q*  and  solve  for  the  matrix 


Nil 
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which  by  our  assumption  is  now  time- independent,  from  equation  (30),  namely 

R*  =  Q*  -  ZP  =  0 

That  describes  the  steady  state  flow  discussed  above.  We  then  derive 
information  about  the  values  of  P(t)  and  P(t)  ::  — at  various 
times  t  from  clinical  data  of  simultaneous  pulse  wave  recordings  at 
the  different  compartments  (such  as  from  Hamit  ct  al.,  loo  oft.). 

Equation  (30)  now  yields  a  n-set  of  relations  for  the  compliances  for  each 
time  value  t  which  is  of  the  form 


(31) 
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Applying  the  Gauss-Markov  theorem,  the  c  values  are  obtained  by  an 
assortment  of  the  set  of  informations  at  all  (n)  given  time  observations, 


namely 


.  W  *  _  7  .  rp 

c  =  nr  iu  ir b 


r‘ £  J  n  x  /> 
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From  the  clinical  measurements  of  Hamit  nt  al.,  (1965),  the  five 
compartmental  compliances  listed  in  equation  (33)  have  been  calculated 


for  the  values 


(X  = 

1000,  6  = 

1/1000. 

They  are 

*T> 

Co 

=  0.0155  ml /mm  Hg 

ironist,  for ■ 

CCF 

=  0.0364 

II 

CBV 

=  0.3180 

II 

CFS 

=  0.0334 

II 

(eonr.t.  for  B  ) 

CFB 

=  0.1830 

II 

At  a  later  stage,  the  numerical  solution  has  been  extended  to  include 

the  general  case  of  C..=C..( a.&).  The  three  dimensional  C..(a,\i) 

ij  ij  ij 

plottings  will  be  reported  separately.  This  concludes  the  calibration 
of  the  compartmental  model  for  the  cerebrovascular  flow  system. 


EVALUATION  OF  THE  COMPARTMENTAL  PULSE  WAVE  FORMS 

Within  the  framework  of  this  work,  which  discusses  the  compartmental 
cerebrovascular  flow  problem  with  constant  resistance  H..  and  compliances 

I'd 

C..  namely,  the  general  linear  problem,  once  the  model  has  been  calibrated. 

I'd 

it  is  possible  to  obtain  solutions  for  the  pressure-time  variations  P(t) 
at  all  the  compartments. 

Solution  of  the  matrix  equation  (12)  in  incremental  form  is  given  by 

p(t+ht)  =  exp( -S'  t)-f 1(A+  (39) 


Here,  denotes  a  time  step  between  the  frontier  time  level 

and  backtime  level  t ^  .  The  matrix  t  equals 


k  =  M-r  7 


(40) 
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The  rational  expansion  of  exp(-K)  leads  to  the  following  formula 

I  -  (1-6)k 

exp(-K)  ~  . . .  (0  <  0  <  1)  (41) 

=  -  I  -  0  k  -- 

where  I  is  the  6*6  unit  matrix.  A  substitution  of  equation  (41) 
into  equation  (40)  now  reads 

-(t+[\t)-  (L  +  -  h-o;k] ir(t)-  '£IQ'\  +  (42) 

The  coefficient  0  controls  the  type  of  solution  which  evolves  in  time. 
When  0  =  0  ,  we  have  an  explicit  scheme;  0  =  1  -  an  implicit  scheme,  and 
0  <  0  <  1  —  a  mixed  scheme.  Thus,  with  the  choice  of  0  ,  equations 
(18),  (21)  or  (24),  (38)  and  (42)  constitute  the  complete  solution  of 
the  non-steady,  compartmental  cerebrovascular  flow  problem  with  constant 
resistances  and  compliances.  In  addition,  being  phrased  as  a  first-order 
differential  matrix  equation  -  Eq .(12)  -  sufficient  boundary  conditions 
ought  to  be  given  to  enable  us  derive  complete  expressions  for  the 
compartmental  pressure  waves.  Unfortunately,  because  of  lack  of  clinical 
data,  rigour  has  to  give  way  to  speculations  from  this  stage  onward. 

We  can  only  express  the  hope  that  the  information  gap  will  close  itself 
before  long. 
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ABSTRACT 

A  lumped  parameter  compar tmentai  model  for  the  non-steady 
flow  of  the  cerebrovascular  fluid  is  constructed.  The  model 
assumes  constant  resistances  that  relate  fluid  flux  to  pressure 
gradients,  and  compliances  between  compartments  that  relate  fluid 
accumulation  to  rate  of  pressure  changes.  Resistences  are 
evaluated  by  using  mean  values  of  artery  and  cerebrospinal  fluid 
(CSF)  fluxes  and  mean  compar tmentai  pressures.  Compliances  are 
then  evaluated  from  clinical  data  of  simultaneous  pulse  wave 
recordings  in  the  different  compartments.  Estimate  of  the 
average  CSF  compartmenial  deformation,  based  on  the  compliance 
between  the  CSF  and  brain  tissue  compartments,  proves  to  be  of 
the  order  of  magnitude  of  actual  experimental  measurements. 


introduction 


The  lumped-parameter  c o mpa r c m en t a  1  model  of  the 
cerebrovascular  system  is  the  first  step  towards  the  construction 
of  a  more  comprehensive  model  of  the  intracranial  fluid  system  . 
The  corapartraental  approach  assumes  that  the  intracranial  content 
may  be  divided  into  a  number  of  units,  or  compartments,  the 
behavior  of  each  of  which  is  represented  by  a  single  value  of 
pressure  and  by  values  of  flux  exchanged  with  adjacent 
compartments.  All  these  values  may  be  time  dependent,  but  they 
don't  vary  in  space.  The  resistance  to  the  flow  from  a 
compartment  to  an  adjacent  one  is  lumped  at  the  boundary 
between  the  two  compartments.  Likewise,  the  integrated  change  in 
volume  of  two  adjacent  compartments  due  to  the  movement  of  their 
common  boundary  in  response  to  a  pressure  difference,  is 
represented  as  a  property  called  compliance  that  is  assigned  to 
the  boundary  between  to  two  adjacent  interacting  compartments. 

Monro's  (1783)  first  model  of  the  intracranial  cavity  was 
b  i -compart mental :  Brain  fluid  and  blood,  as  two  almost- 
incompressible  fluid  phases.  Kellie  (1824)  modified  this  model 
by  assuming  three,  instead  of  two,  material  compartments; 
arteries,  veins  and  brain  tissue.  The  Monro-Kellie  models 
prevailed  to  this  century  and  were  modified  in  stages  only  in 
recent,  years.  The  number  of  compartments  was  increased  to  six: 
arteries,  capillaries,  venous,  venous  sinus,  jugular  bulb  and 
cerebrospinal  (CSF)  (Agarwal,  1971). 

Yet,  in  all  these  models,  the  fluid  itself  remained 
incompressible.  More  recent  approaches  relaxed  the  latter 
constraint  and  the  fluid  was  taken  as  (linearly)  compressible. 
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As  to  the  brain  tissue,  it  is  essentially  a  raultiphasic  material 
(e.g.,  brain  tissue,  blood  and  C3F)  continuum.  Experimental 
results  show  that  its  behavior  is  inelastic  or,  alternatively , 
that  the  corresponding  compliance  is  non-linear.  To  overcome 
the  "non-linearity”  of  a  single  coefficient,  the  tissue  is 
assumed  to  be  a  rather  complex,  single-phase,  multi-parameter 
viscoelastic  material  e.g.,  one,  whose  constitutive  relation 
involves  four  viscoelastic  coefficients.  (Panidi  and  Advam 
19 7 8 ) .  In  the  reports  of  these  investigations,  although  not 
explicitly  stated,  the  model  of  the  intracranial  content 
returned  to  be  b i -c ompar t men t a 1 :  The  CSF  compartment  and  "all 
the  rest",  or  the  vascular  compartment  and  the  rest,  etc.  In 
most  cases,  even  for  the  multicompartmental  model  of  the 
cerebrovascular  fluid  system  (e.g.,  Agarwal,  1971),  no  numerical 
calculations  were  presented  and  the  exposition  of  the  subject 
remained  theoretical. 

Our  first  objective,  therefore,  was  to  develop  an  N- 
compartmental  model  that  can  yield  numerical  values  of  the 
various  state  variables  (e.g.,  pressure).  So  far,  we  have 
successfully  achieved  (Sorek  et  al.  1986a)  this  objective  for 
the  general  linear1  problem,  assuming  constant  compliances  and 
resistances.  Physiologically,  the  assumption  of  linearity 
corresponds  to  passive  states  in  which  the  sensory  and 
endocrinological  biocontrol  mechanisms  have  but  little  effect  on 
the  resistances  and  the  compliances.  Basically,  homeostatis  is 
non-linear  by  virtue  of  the  feedback  mechanism  which  senses  the 
deviations  in  the  values  of  resistances  and  compliances  and  acts 
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to  restore  them  to  normal  values.  Nevertheless,  one  must  first 
construct  and  solve  a  linear  model  and  employ  perturbation 
techniques  in  order  to  derive  solutions  for  the  generalized  non¬ 
linear  problem.  Accordingly,  we  have  first  constructed  and 
solved  a  linear  model  that  described  non-steady  flow  in  a 
lumped-parameter  N-compar tmental  model  of  the  cerebrovascular 
fluid  system  with  constant  resistances  and  compliances. 

The  compartmental  model  involves  a  number  of  resistances  and 
compliances,  the  values  of  which  must  be  known  before  the  model 
can  be  employed  in  predicting  pressure  and  flux  changes.  In  the 
first  paper  (Sorek  et  al.  1966a)  it  was  indeed  assumed  that  the 
values  compliance  and  resistances  including,  as  we  shall  see  also 
here,  certain  ratios  between  pairs  of  resistances  were  known. 
The  objective  of  the  present  paper  is  to  present  a  methodology 
for  estimating  the  values  of  the  variouis  compliances  and 
resistances  and  the  above  mentioned  ratios. 

1.  THE  COMPARTMENTAL  BALANCE  EQUATION 

The  governing  equation  for  the  lumped-parameter 
compartmental  model  describe  the  balance  of  mass  and  the  balance 
of  linear  momentum  for  each  compartment.  Essentially,  each  such 
equation  states  that  the  temporal  rate  of  change  of  either  the 
fluid  mass,  or  its  momentum  in  a  compartment,  is  equal  to  the 
amount  of  net  influx  of  that  quantity  through  the  compartmental 
boundaries  plus  the  external  sources.  The  mass  balance  of  the  n- 
th  compartment,  surrounded  by  a  numbr  of  compartments  denoted  by 
m=l,2,...,  can,  therefore,  be  written  in  the  form 


dV 

— ~  =  T  q  +  Q 
dt  (£}  mn  n 
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(1) 


where  q^  ^(-q^  )  denotes  the  flux  from  compartment  n  to  m, 

denotes  external  sources  in  the  n-th  compartment  and  V  is  its 

n 

volume . 

The  flux  q  can  be  expressed  in  ter- ms  of  the  difference 

nm 

in  pressure,  P  (=P  -P  )  between  the  n  and  the  m  compartment,  and 
nm  n  m 

a  conductance  Z  (reciprocal  of  the  resistance  H  ),  in  the 

nm  nm 

f  orm 


q 


nm 


P 

nm 

R 

nm 


Z  P 
nm  nm 


(2) 


The  change  in  volume,  a V  is  produced  by  the  change  in 

n 

the  pressure  differences,  AP  ,  in  adjacent  compartments,  taking 

nm 

into  account  the  presence  of  compliances,  C  ,  between  these 

nm 

cells 


« 


ye  ap 

k  nm 


(m) 


nm  nm 


(3) 


Together,  we  obtain  for  the  n-th  compartment,  a  mass  balance 
of  the  form 

dP 


l  c 

(m) 


nm 


nm  dt 


I  z  p  =  q 

,  nm  nm  n 

(m) 


(4) 


Another  compact  form  of  this  equation  for  ail  cells 
simultaneously  is 


dP 

SdT 


i 


♦  z  p 


9 


(5) 


62 


where  £(t)  is  the  time-dependent  (Nxl)  pressure  vector;  Q(t)  the 
source  (Nxl)  flux  vector,  and  Z  and  £  are  the  (NxN)  conductance 
and  the  compliance  matrices,  respectively. 

2.  PARAMETER  ESTIMATION 

The  corapartmental  balance  equations  for  fluid  mass,  written 
in  the  compact  form  of  (5),  involve  conductivities  and 
compliances,  expressed  by  the  matrices  Z  and  C^.  In  the  present 
work,  these  parameters  are  assumed  to  be  constant. 

To  predict  the  pressure  (and  flux)  response  of  the  model  to 
external  changes,  the  values  of  the  parameters  C  and  Z  must  be 
known.  In  order  to  estimate  them,  we  need  measured  values  of 
pressure  in  the  various  compartments  at  a  sufficient  number  of 
points  in  time. 

Because  C  and  Z  are  constant,  and  p(t)  is  cyclic,  by  taking 
a  temporal  average  of  (5),  (i.e.,  integration  over  a  period  of 

time  divided  by  the  period),  we  obtain, 

-*  *  . , . 

P  Z  =  Q  (o) 

• 

where  (  )  denotes  a  temporal  rate  of  change,  (  )  denotes  the 
difference  between  adjacent  compartments  and  p*  and  Q*  denote 
mean  values  of  pressure  differences  and  source  fluxes, 
respectively.  We  note  that  (6)  is  a  quasi-steady  state  equation. 
With  known  values  of  P*  and  Q*,  we  solve  (6)  for  Z. 

Given  further _ information  from  clinical  data  e.g.,  Hamit  et 
al.  196b),  of  simultaneous  pulse  wave  recordings  P(t)  and  P(t) 
(idP(t)/dt),  at  the  different  compartments  and  at  various  times, 
t  ,  k  =  l,2,...K,  equations  (3)  and  (4)  now  yield  a  K-set  of 


relations  for  the  compliances. 


rc  =  q  -  zr  ( 7 ) 

By  the  Gauss-Markov  theorem,  the  C  values  are  derived  as  an 
assortment  of  the  set  of  informations  at  all  K  time  observations 
(Gorek  et  al.,  1986a). 

This  concludes,  at  least  formally,  the  (inverse)  process  for 
identifying  parameters  appearing  in  equation  (5)  for  a  N- 
compartment  model. 

3.  EVALUATION  OF  MODEL  REGIGTArtC&G  AnD  COMPLIANCES 

Let  us  determine  the  values  of  C  and  Z  in  the  case  of  a 
seven-compartment  model,  namely  for  N=7. 

Figure  1  shows  the  model  consisting  of  the  following  compartments: 

Fig.  1 

arterial  (A),  capillary  (C),  cerebrospinal  fluid  (F),  brain 
tissue  (B),  venous  (V),  venous  sinus  (S)  and  the  jugular  bulb 
(J).  The  (lumped)  resistances  are:  between  the  artery  and 

capillary  compartments  (R^C);  the  capillary  and  cerebrospinal 
fluid  compartments  ( Rcf  > >  the  capillary  and  brain  tissue 
compartments  (R^b),  the  capillary  and  vein  compartments  (R  )» 
the  brain  tissue  and  vein  compartments  (Rgv  )»  cerebrospinal 

fluid  and  brain  compartments  (RpB),  the  vein  and  venous  sit. us 
compartments  (R  )>  the  cerebrospinal  fluid  and  the  venous  sinus 

compartments  (Rps)f  between  the  venous  sinus  and  the  jugular  bulb 
compartments  ( Rv ,  ) 5  altogether  nine  resistances.  In  Figure  1, 
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4 

t 


the  capillary  compartment,  is  divided  into:  the  choroid  plexuses 
-  those  tufts  of  small  capillary  vessels  inside  each  of  the  four 
ventricles  -and  the  capillary  system  outside  the  ventricles. 
However,  in  the  equations,  only  the  combination  in  parallel  of 
che  conductances  z;c  ♦  z"c  -  zAC  and  z^  .  z“v  -  zcv  appear,  s0 

that  only  the  combined  resistances  RAC  and  F^v  ~  into  and  out 
of  the  capillaries  -  are  included  in  the  model. 

The  resistances  ,  RCf.  and  RfB  are  identified  as  the 
lumped  blood-brain  barrier;  the  lumped  blood-cerebrospinal  fluid 
barrier,  and  the  lumped  cerebrospinal  fluid  -  brain  barrier, 
respectively  . 

We  recall  that  the  compliance  elements,  Cnm>  indicate  that 
an  increase  in  volume  of  one  compartment  equals  the  volume  of  the 
"cup"  formed  by  the  deformed  membrane.  This  volume,  in  turn, 
equals  the  volume  displaced  from  the  neighboring  compartments, 
all  this  within  the  rigid  container  of  the  skull  bones  (the 
Monro-Kellie  doctrine). 

In  the  non-steady  state,  which  takes  into  account  the 
de f or mab i 1 i 1 ty  of  the  compartments,  we  first  introduce  a 

compliance  element  between  the  artery  and  the  brain  tissue 
compartments  it  represents  the  overall  pulsatory  effect  of  the 
arteries  on  the  brain  tissue.  Next,  the  capillary  system  is 
considered'  non-def ormable,  so  that  no  compliance  is  introductea 
between  this  compartment  and  any  of  its  neighbors.  The  choroid 
plexuses,  however,  although  capillary  in  nature,  possess  other 
material  properties.  Hence,  they  can  ,  and  in  fact  do,  convey 
pulsations  to  the  CSF  system  (Bering,  1955).  Accordingly,  a 


introduced  between  them.  Furthermore,  the 
share  common  boundaries  -  at  the 
subarachnoidal  space  -  which  are 
element  Cft  Is,  therefore,  inserted 
Finally,  to  account  for  observed  sharp  drop  in 
the  cardiovascular  passage,  additional  compliances 
are  inserted  between  the  brain  tissue  and  venous 
and  between  the  CSF  and  venous  sinus  compartments, 
Altogether,  in  our  presentation,  we  assume  five 
adjacent  elements  of  the  cerebrovascular 


compl iance  C  is 
CSF  system  and  the  brain  tissue 
ventricles  and  along  the 
deformable.  A  compliance 
between  the  two. 
pressure  along 

<-'uv  an<*  Cps 
compartments 
respectively, 
compliances  between 
fluid  system. 

The  mechanical  properties 
are  symmetric  with  respect  to  the 
compartment  and  its  neighbor,  i.e, 


R 

=  R  , 

R  = 

ft 

AC 

CA 

CF 

FC 

C 

=  c  , 

C  = 

C 

AB 

BA  ’ 

CF 

FC 

of  resistances  and  compliances 
change  of  direction  between  one 
in  formulae 
etc  . 
etc . 


All  H  ’s  and  C  's  are  positive, 
am  am 

We  adopt  the  value  of  mean  pressures  along  the 
cerbrovascular  fluid  system  according  to  the  pressure  profile  of 
the  car d i o vascu 1 ar  system  cited  in  the  literature  (e.g.  in 
Guyton,  Ch.  14,  1969).  In  the  arteries,  the  average  pressure  - 
between  systole  and  diastole  -  is  P*  =100mmHg.  It  drops  to 

A 

P*  =30mmHg  in  the  capillaries,  including  the  choroid  plexuses; 

P*  =10mmHg  in  the  CSF  system,  P*  =9.5mmHg  in  the  brain  tissue, 
F  B 

P*  =9mmHg  in  the  venous  system,  P*  =8mmHg  in  the  sinuses  and 
P*  =2mmHg  in  the  larger  juglar  veins  and  in  the  spinal  leading 
into  the  vena  cava  (Fig.  1).  The  mean  values  of  the  injected  and 
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ejected  fluxs  at  the  artery  and  jugular  bulb  are 
Q*A=Q*j  =750ral/min. 

Altogether,  nine  values  have  to  be  determined.  However,  the 
matrix  equation  (6)  comprises  only  six  independent  balance 
equations  for  the  various  compartments  of  the  cerbrovascular 
model.  Thus,  the  redundancy  of  the  system  is  three,  and  three 
additional  conditions  are  needed  to  solve  the  set.  One  of  them 
is  the  mean  flux  (from  choriod  plexus  to  the  CSF  ventricles) 

Q*F=0.3ml/min,  which  can  be  taken  as  pivot  value  with  high 
credibility  from  the  literature,  (Cutaler,  1968). 

Thus,  the  Zqp  value  can  be  evaluated  from  the  expression 


* 


* 

z  p 

CF  CF 


(8) 


Two  more  conditions  have  to  be  stipulated  for  resistences, 
based  on  the  existing  physiological  data.  We  now  introduce  the 


scalar  coefficients. 


rfb  _  !vs 

RVS  ZFB 


(9) 


RFB  =  fcB 

rcb  zfb 


(10) 


Here,  a  indicates  the  ratio  of  the  resistance  of  the 
cerebrospinal  fluid-brain  barrier  to  the  vein-venous  sinus 
resistance,  whereas  B  is  the  ratio  of  resistance  of  the 
cerebrospinal  fluid-brain  barrier’  to  that  of  the  blood-brain 
barrier.  Thus,  equations  (6)  to  (10)  allow  a  complete  solution 
for  the  resistences  R  (a,  B  )  and  compliances  Cnm  (a>B)  with  the 
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values  of  a  and  6.  Figures  2  and  3  describe  an  example  of  the 
surfaces  H  (  a,  6)  and  C 

L  V  BV 

respectively 

Figs.  2,3 

The  figures  demonstrate  zones  of  a  and  g  that  generate 
unacceptable  values  (e.g.  negative  Z’s  as  in  Fig.  2)  of 
resistances  and  compliances  and  zones  of  high  sensitivity  of  the 
resulting  Z's  and  C's  to  small  changes  of  a  and  8.  Our  choice 
is,  therefore,  to  rely  on  <*  and  8  values  that  generate  stable 
Z's  and  C's. 

There  are  almost  no  data  about  the  physiological,  or 
pathological,  ranges  of  a  and  8.  It  was  found  that  when  a.  and  8 
are  in  the  range  of  1/1,000  and  10,000,  respectively  the 
resistances  and  compliances  meet  the  desired  aforementioned 
criteria.  Hence  for  6  =10  and  <*=10  ,  equations  (4)  to  (8) 

result  in  the  following  values 

^AC=0 .0933mmHg/ml/min .  Hvg=0.00ll  mmHg/ml/min. 

ft =  66  .667 mmHg/ml/min  .  H  =  7  .6187 mmHg/ml/min  . 

RCV  =  0  .028ramHg/ml/min .  Rgj  =0  .0080mmHg/ml/min . 

=  13338  .OmmHg/ml/min  .  ^FB  =  ^  •338mmHg/ml/min  * 

HBV  =  3  . 33romHg/ml/min . 

and 

^aB  =  0  .0012ml/mmHg  C  =0  .0494ml/mmHg 

r  b 

C,  =0  .0357ml/mmHg  C  =C- .2093wl/mmHg 

C  r  FB 

'-'BVsO  .3746ml/mmHg 
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With  the  above  value  of  C  t->  a  >  we  can  nQw  assess  the  average 
deformatibn  of  the  C3F  compartment.  Let  us  assume  a  spherical 
conf iguration  of  this  compartment,  with  a  mean  diameter  r^,  xts 

/  n 

volume,  V,,,  and  surface  area,  are  given  by  v„  =  -rnr', 

r  r  r  i  r 

anJ  s(,  -  .  dy  virtue  of  equation  (2),  we  may  thus  express 

the  change  in  V  ,  by 
t' 


AVF=  VrF 


=  C  F 
FB  FB 


Ul) 


According  to  Hakim  et  al.  (1976),  the  mean  diameter  of  the 
CL F  compartment  is  ...  -  }  *  Thus,  in  view  of  the  mean 

h  4  2  if 

•k 

pressure  difference  i’K1J  =  0.bmmllg,  the  compliance  value 
=0.C093ml/mmHg,  and  from  equation  (11),  we  obtain 


*'F=  16  .07*  10~b  mm. 


This  estimate  of  displacement  of  the  CLF  compartment 
boundaries  is  consistent  with  measurements  done  by  Allen  et  al. 
(  1983)  . 

Finally,  we  wish  to  emphasize  that  the  model  approach 
p resented  here  (see  also  Sorek  et  al.,  1986b,  Kami  et  al., 
I9du  )  constitutes  a  methodology  that  mn  be  implemented  to 
various  compartmental  schemes  representing  different  aspects  of 
cl i n i cal  data . 
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PUMP  OUTLET 


Lumped-parameter,  seven-compartmental  model  of  the  intracranial 
cerebrovascular  fluid  system  ( R-resistance ,  C- compliance; 

C  3  pressure  (mm  Hg);  (  )  flow  (ml /min) ;  <  >  volume  (ml)). 


RESISTANCE  R(CV)  -  VS  -  BETA  Sc  ALPHA 

f  XT  T'  n  n  n  a  rm  r-*  * 

/WJBf  TA=2100  IN  ALP  HA  =  4000 


INCREMENTS: 


Fig.  2:  SURFACE  OF  R  (<*,£) 
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INTRACRANIAL  COMPARTMENTAL  PULSE  WAVE  SIMULATION 
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Department  of  Biomedical  Engineering  and  Department  of  Civil 
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Abstract 

The  general  solution  of  the  linear  compartmental  model  for 
the  cerebrovascular  fluid  system  with  constant  resistances  and 
compliances,  predicts  the  pressure  waves  in  the  compartments  in 
response  to  an  input  pulse  wave-arterial  and/or  jugular.  Results 
are  shown  for  a  seven-compartment  model  and  for  a  sinusoidal 
arterial  pulse  wave  at  a  frequency  of  1  Hz,  with  and  without 
fluid  drainage  from  the  CSF-compartment . 


Keywords:  Brain  tissue,  Compartmental  model,  Intracranial 

pressure,  Pulse  wave. 
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INTRODUCTION 


We  chose  to  simulate  the  fluid  regime  in  the  brain  by  a 
seven-compartment  model  (n  =  7 )  shown  in  Figure  1.  It  is  based  on 
Agarwal's  (1971)  six-compartment  model:  artery  (A),  capillary 
(C),  cerebrospinal  fluid  (F),  vein  (V),  venous  sinus  (S)  and 
jugular  bulb  (J),  to  which  we  added  the  brain  tissue  compartment 
(B).  The  model  assumes  nine  resistances  and  five  compliances  as 
shown  (Sorek  et  al,  1986a  and  b)  in  the  figure  (note  that  the 
capillary  compartment  is  divided  into  two  parallel  one,  the 
choroid  plexuses  inside  the  ceberal  ventricles  and  the  rest  of 
the  capillary  system  outside  the  ventricles.  Thus,  the  combined 
resistances  R^q  and  r  enter  into  the  calculations). 

L  V 

The  mean  compar tmental  pressures  are  also  shown  on  Fig.  1 

p*  (average  between  systole  and  diastole)  =100mmHg,  P*  =30mmHg, 
A  u 

P*  =  lOmmHg,  P£  =9.5mmHg,  P$  =9 mmHg,  P*  =8mmHg  and  P*  =2mmHg. 
The  mean  arterial  flux  is  Q*  =750ml/min. 

The  intracranial,  compart  mental  fluid  flow  of  the 
cerebravascular  system,  is  governed  by  the  fluid  balance  equation 
(Sorek  et  al  .  1986a)  . 


dp 

£h7  =  9 


(1) 


This  is  an  inhomogeneous,  ordinary  differential  matrix 
equation  of  the  first  rank  with  respect  to  time  (t).  Here,  C  is 
the  NxN  matrix  of  the  intercompartmental  compliances;  Z  is  the 
NxN  matrix  of  the  reciprocal  intercompartmental  resistances,  P 
is  the  N x 1  matrix  of  the  compartmental  pressures,  Q-  the  Nxl 
matrix  of  the  compartmental  sources  which  add  fluids  to  the 
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intracranial  domain  from  outside.  The  fluid  density  is  assumed 
constant . 

Equation  (1)  represents  here  six  linear  balance  equations, 
one  for  each  compartment,  except  the  Jugular  bulb,  in  the  six 
variables  (Pa  p  p  p  p  p  \ 

"  ’  B  ’  C’  F  ’  *V  »  *S  >’ 

The  components  of  the  Q  matrix  are:  QA>  0,0,0,0,,-q  where 
Qa  is  the  flux  input  to  the  artery  and  Q(  Tlux  outflow 

from  the  Jugular  bulb. 

Equation  (1)  is  subject  to  initial  pressure  values 


p  =  p 

==(  t  =0)  -0 


( a) 


Basically,  the  problem  is  non-linear.  By  virtue  of  the 
homeostatic  biocontrol  (sensory  and  endocrinological)  feedback 
mechanisms,  the  coefficient  matrices  £  and  Z  can,  and  in  fact  are 
functions  cf  the  pressure  vector  £  and  of  time,  t.  There  are, 
however,  some  passive  cases  which  can  be  approximated  by  the 
linear  problem,  namely  when  C  and  Z  assume  constant  values. 
Moreover,  even  if  we  deal  with  the  non-linear  problem,  the 
solution  passes  through  the  linear  case.  In  what  follows  we  will 
concentrate  on  the  solution  of  the  linear-  problem,  i.e.,  with 
constant  C  and  Z . 

In  what  follows,  we  assume  that  the  actual  values  of  the  Z 
and  C  components  are  known.  These  are  obtained  by  solving  the 
appropriate  parameter  estimation,  or  inverse,  problem  (Sorek  et 

al .,  1986b)  . 

Actually,  in  view  of  the  Monro-Kellie  doctrine  which  states 


that  the  sum  of  all  compartmental  volume  changes  equals  zero, 
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N 

l  AV  =  0  v3) 

n=l  n 

Only  five  independent  algebraic  equations  are  represented  by 

(1)  . 

In  other  words,  if  one  compartment  changes  volume,  it  has  to 
be  at  the  expense  of  the  volume  of  its  neighbors,  all  within  the 
cranial  vault  assumed  to  be  absolutely  rigid.  This 
intercompartmental  relation  is  reflected  in  the  expressions  for  Z 
and  C  . 


I  C,,  ■  Sc,,  -  0 


1J 


1J 


m 


l  z. .  =  y  z. .  =  o 

^11  i  i 


ij 


j  1 

This  means  that  we  need  external  information  on  the  pressure 
in  one  of  the  compartments.  Because  the  pressure  artery  can  be 
measured,  (perhaps  in  the  future  by  non-invasive  methods),  we 
choose  to  specify  it  as  the  additional  information. 

In  the  following  sections,  we  bring  the  numerical  results 
for  the  compartmental  pressure  waves  that  result  from  an  input 
arterial  sine  wave  with  and  without  a  sink  (drainage)  in  the 
Cerebrospinal  Fluid  (F)  compartment.  We  will  also  discuss  the 
physiological  and  clinical  significance  of  some  of  these  results. 


2 .  COMPARTMENTAL  PULSE  WAVE  FORMS 

The  choice  of  a  sine  wave  as  the  input  function  QA  was  macie 
for  two  reasons.  First,  some  of  the  observed  pressure  wave 
forms  are  indeed  close  to  the  sine  waves.  Secondly,  if  they  do 
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not  correspond  to  sine  waves,  but  are  still  periodical,  then  by 
Fourier's  analysis  they  can  be  described  as  linear  combinations 
of  sine  waves  of  different  frequency,  but  of  the  same 
periodicity.  We  also  choose  the  frequency  of  1  Hz  (1  cycle  per 
second),  as  base  frequency  for  the  calculations.  This  frequency 
is  close  to  the  cardiac  frequency.  In  the  classification  of  the 
intracranial  pressure  wave  forms,  it  corresponds  to  the  Lundberg 
A-wave  (1974),  also  termed  "pulse  wave". 

Pa  =  P,  ♦  P.  sin  (wt) 

A  A  A 

Where  P*  =100mmHg  is  the  average  arterial  pressure, 

I 

PA=20mmHg  is  the  amplitude  of  the  arterial  pulse  wave 
between  systole  and  diastole,  and u  = £* Hz ,  namely  the  period  of 
the  wave  is  one  second  upper  part  of  Fig.  p_ . 

Based  on  the  above  considerations  on  l\(t)  >  the  arterial 

A 

fluctuating  flux  was  chosen 


Qa  =  QA  +  Qa  sin  (wt )  (6 ) 

where  Q*  =750ml/min,  and  q'  =125ml/min  for  U-uh-h  and 

A  A 

®A  =100ml/min  for  ir<wt<2ii 

Jagular  bulb  involves  no  compliance  (i.e.,  its  walls  are 
practically  non-de  f  ormab  le)  .  Hence  the  fluid  mass,  balance 
equation  for  it  is 


VW  - 


(7) 
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Due  to  the  absence  of  compliances  in  the  Jagular  buib,  we 
assume  that  the  pressure  in  it  under  unsteady  conditions  (i.e., 
in  reality)  is  the  same  as  under  quasi  steady  conditions.  The 
latter  are  defined  by 


Z  P 


(t) 


{(t) 


(8) 


i.e.,  deleting  the  effect  of  cell  compliances.  Under  quasi 
steady  state  conditions 

QA<t>  =  Qj(t)  ('.)) 

We  take  advantage  of  this  conclusion  by  combining  (7) 
through  (9)  to  yield  the  pressure,  P;(l)>  in  the  jaguiar.  bulb 

compartment.  Furthermore,  knowing  the  pressure,  P  »  we  employ 
(7)  to  express  Qj  in  terms  of  Pj  and  substitute  this  expression 
into  Qj(0  given  in  (1)  and  solve  for  the  pressure  variations 
P(t)  in  the  non-steady  case. 

We  should  emphasize  that  (9)  is  valid  only  for  a  quasi¬ 
steady  state.  Under  unsteady  state  conditions,  QA< 1  )=  q  (t). 

Once  a  solution  for  Ps(t),  among  other  pressure  values,  is 
obtained  one  may  insert  it  into  (7),  which  is  valid  also  for 

unsteady  state,  to  obtain  Qj(t)  under  unsteady  flow  conditions. 
It  is  also  possible  to  calculate  the  arterial  flux  by  (eqs.  (l),  (■>)) 

*  •  dPB 

QA  ’  VVV,i1Vw“)*c».“"l"t>l*t„ir  (10) 


obtained  from  (1),  that  also  takes  into  account  the  compliances 
in  all  of  the  compartments. 
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3 .  EXAMPLES 

A  number  of  computer  solutions  were  derived  for  cases  of 
interest.  In  all  cases,  the  following  values  of  resistances  and 
compliances  were  employed  (Sorek  et  al .  1986b). 


R  =0  .0933mmHg/ml/min . 

AC 

u  =0  .0013mmHg/ml/min . 

VS 

rcf  =66  .667mmHg/ml/min  . 

R  =7 .6 l37mmHg/ml/mi n . 

R(  V  =0 .028mmHg/ml/min  . 

, =0 .0080 mmHg /ml /mi n  . 

R  ...  =13338  .OmmHg/ml/min  . 

LB 

R  =13 .338mmHg/ml/min . 

Rbv=3-33  mmHg/ml/min 

and 

t:  =0 .0012ml/mmHg 

Ail 

(  (  ^  =0  ,0494ml/mmHg 

c  =0  -0357ml /mmHg 

a 

C =0  .02093ml /mmHg 

c  =0 ,3746ml/mmHg 

EXAMPLE  1; 

Given  the  phasic  pressure  PA(t),  as  expressed  by  eqn. 
(5),  and  the  arterial  flux,  Qa)  >  expressed  by  eqn.  (6),  the 

A 

Jugular  phasic  pressure,  Pj(t)  and  the  corresponding  pressures 
P  (t),  P  (t),  P  a)  ,  P„(t)  and  P  (t)  were  determined.  The 

C.  d  r  o  v 

results  for  Pg(t)  are  shown  in  Fig.  2 .  l'he  other  pressures  are 

shown  in  Fig.  3- 

It  is  of  interest  to  note  (in  Fig.  3)  the  negative  values 
0g  p  (t)  during  part  of  the  cycle,  as  a  result  of  the  suction 
action  of  the  cardiac  system.  This  is  consistent  with  clinical 


observations . 
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Figure  4  shows  the  calculated  Jugular  outflow  and  arterial 
inflow.  As  may  be  seen  from  the  figure,  the  two  fluxes  are 
equal  to  each  other.  The  time  integrated  difference  over  half  a 
cycle  indicates  no  fluid  quantity  stored  and  then  released  from 
storage  during  each  cycle. 

EXAMPLE  Z: 

By  introducing  a  shunt  into  the  CSF  ventricles  (F 
compartment),  CSF  can  be  drained  out.  We  impose  drainage  rates 

of  .2,  .8,  1.4ml/min.  Figure  5  shows  the  resulting  pressure, 

Pp(t)  .  As  can  be  seen,  an  excessive  drainage  will  lead  to  a 

strong  distortion  of  the  sinusoidal  wave  configuration. 

EXAMPLE  3: 

Upon  imposing  a  constant  pressure  on  the  brain  tissue, 

P  =  9 .5 m mHg ,  Figure  6  shows  that  the  pressure  curve  in  the  F- 

B 

compartment  will  be  highly  decayed.  This  fact  emphasis  the 
importance  of  the  pulsational  behavior  of  the  pressure  within 
each  compartment,  of  the  cerebral  system. 

SUMMARY 

The  interaction  of  the  various  components  of  the 
cerebrovascular  fluid  system  is  representable  by  a  lumped- 
parameter  compar tmental  model  which  takes  into  account  the 
pressure  and  volume  changes  between  the  compartments. 

Presently,  only  the  linear  problem  is  discussed  for  which 
the  resistances  and  compliances  are  assumed  constant.  The 
general  solution  of  the  linear  problem  enables  us  to  derive  the 
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non  steady  pressure  waves  in  all  the  compartments,  also  with  sink 
(drainage)  function  systems. 

The  computer  has  been  programmed  for  the  general  case  of  an 
N-compartmental  model.  The  results  brought  here  are  for  a  seven- 
compar tmental  (N=7)  model:  arterial,  capillary  (including  the 
choroid  plexuses),  cerebrospinal  fluid,  brain-tissue,  venous- 
sinus  and  jugular  bulb. 

The  simulation  of  the  compartmental  pressure  waves  is  for  an 
input  sinusoidal  arterial  pressure  wave,  also  termed  "pulse 
wave",  of  1  Hz  frequency. 

« 

The  linear  assumptions  correspond  fairly  well  to  some 
passive  states  of  long-standing  coma  and  chronic  neurologic 
cases.  Even  when  the  problem  becomes  non-linear  and  the  feedback 
effects  of  the  bio-control  mechanisms  -  sensory  and 
endocrinological  -  have  to  be  incorporated  into  the  model,  the 
linear  solution  will  again  be  the  first  stage  upon  which  the 
nonlinear  effects  are  then  superimposed.  We  thus  have  at  our 
disposal,  at  least  as  a  first  step  a  method  of  "computerized 
simulation"  for  the  intracanial  fluid  dynamic  system  the  use  of 
which,  particularly  in  cerebral  intensive  care  unit,  is  self- 
evident  : 
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CAN  N.P.H.  BE  CAUSED  BY  CEREBRAL  SMALL  VESSEL  DISEASE? 

A  new  look  based  on  mathematical  model 

S.  Sorek  D.Sc.* ,  M.  Feinsod  M.D.**  and  J.  Bear  Ph.D.*** 

ABSTRACT 

A  novel  mathematical  model  describing  the  intracranial 
contents  as  lumped  interacting  compartments  is  presented.  The 
model  predicts  pressures  and  fluxes  as  function  of  time  in  the 
various  compartments.  Compart ment al  resistances  and  compliances 
are  evaluated  as  step  functions  of  mean  pressures  and  fluxes 
values.  According  to  this  model,  normal  pressure  hydrocephalus 
may  be  the  result  of  small  vessel  disease  that  abolishes  the 
pressure  gradient  between  the  capillaries  and  brain  tissue. 
Lowering  the  CSF  pressure  as  by  shunting,  restores  the  required 
compart  ment  al  interaction  with  new  values  for  the  resistances  and 
compliances . 

Keywords;  normal  pressure  hydrocephalus,  ventrical-peritoneal 
shunt , mat hemat i cal  model,  intracranial  pressure,  compliances, 
resistances,  fluxes,  CSF  physiology,  capillaries. 
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INTRODUCTION 

The  clinical  features  of  Norma . 
are  well  known  (Fisher,  1977).  The  ir p.'h i 
ventricles  enlarge  in  adult  patients  *  ,-\r.  ,  .. 

intracranial  pressure  is  still  speoulai  ,vr  . 

1977).  Various  tests  were  suggested  in  omie:  .. 

factory  prediction  which  patient  will  here: 
procedure  but  no  one  gained  yet  the  expected  rec»ipr. 

In  this  paper  we  will  report  a  patient  *  „  t  i  , ,  • 
clinical  course  could  be  explained  by  a  p  p  _ : : a : 
mathematical  model  of  inter  relations  between  brain  ct ii.pc.i- 
(Sorek  et  al.  1986a). 

This  model  may  shed  light  on  the  still  obscure  et:n.<.t; 

NPH  . 
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CASE  REPORT 

A.K.  is  a  65  year  old  merchant  started  exhibiting  signs  of 
memory  loss  and  impaired  judgement.  His  appearance  continued  to 
be  immaculate  and  small  talk  did  not  reveal  his  deficit.  CT  Scan 
demonstrated  enlarged  ventricular  system.  The  fourth  ventricle 
seemed  less  involved  than  the  others  CSF  pressure  on  lumbar 
puncture  was  120mmH2<“)  RHISA  cisternography  showed  rapid  entrance 
of  the  isotope  into  the  ventricles,  it  cleared  only  after  more 
than  48  hours. 

A  shunting  procedure  was  suggested  but  the  family  elected  to 
wait.  His  condition  slowly  deteriorated,  his  dementia  became 
overt  and  he  was  confined  to  home.  Six  months  later  another 
consultation  was  seeked  now  because  of  progressive  ataxia. 
Repeat  CT  revealed  further  enlargement  of  the  hydrocephalus. 
Only  4  months  later,  when  the  patient  was  confined  to  bed  due  to 
severe  ataxia,  incontinence  and  speechlessness  did  the  guardian 

permit  operation.  CSF  pressure  at  that  time  was  115mmH2°" 

A  ventruclo-peritoneal  shunt  with  an  opening  pressure  of 

9<JmmH2°  w as  installed.  The  postoperative  course  was  remarkable 
for  the  rapid  return  of  speech,  memory,  ambulation  and 
continence.  A  month  after  the  operation  the  patient  returned  to 
his  business  and  several  weeks  later  reported  success  in 
complicated  financial  considerations  and  decisions.  CT  scan 
demonstrated  small,  well  drained  ventricles. 

Eight  months  later  he  started  to  deteriorate  and  within  3 
weeks  he  was  approaching  his  pre-operative  condition.  CT  scan 
revealed,  again,  enlarged  ventricles  (Fig.  1).  Shunt  malfunction 
was  presumed  but  surgical  revision  failed  to  demonstrate  any 
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obstruction  and  did  not  alter  his  condition.  Only  after 
installment  of  a  new  shunt  system  with  a  low  opening  pressure 

(bOmmH^o)  did  the  patient  improve.  Once  again,  his  recuperation 
was  remarkable;  he  is  back  in  finances  and  doing  well.  CT  scan 

(Fig.  2)  is  evident  for  well  drained  ventricular  system. 


ICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS -1963- a 


DISCUSSION 


In  a  series  of  recent  papers  (Kami  et  al.  1986,  Sorek  et 
al.  1986a, b,c)  a  model  depicting  the  brain  as  an  assembly  of 
interacting  compartments  was  put  forward. 

The  model  (Fig.  3)  solves  the  distribution  of  pressures, 
fluxes .resistances  and  compliances  within  these  compartments. 

In  this  model  we  attempt  to  define  each  intracranial 
structure  as  a  cell.  Seven  such  cells  -  arteries  (A), 
capillaries  and  choroid  plexus  (C),  veins  (V),  venous  sinuses, 
(S),  ventricular  cerebrospinal  fluid  (F),  jugular  bulb  (J)  and 
brain  tissue  (B)  are  lumped  together  and  their  interactions  are 
described  by  a  series  of  flux  balance  equations  (Sorek  et  al. 
1985a)  . 


As  an  example  let  us  consider  the  equation  describing  flux 


balance  for  the  capillary  compartment. 


PA~PC 


P„-P„  P„-P 


P„-P« 


C  F  .  ‘C  ‘B  *C  *V  c  J2_  (P  _p  ) 

RCV  CF  dt  C  F 


where 

PA=  arterial  pressure 
p  =  capillary  pressure 

PF=  CSF  pressure  (ventricular) 

p  =  brain  tissue  pressure 

B 

PV=  venous  pressure 

PAC=  resistance  to  flow  between  arteries  and  capillaries 

p  =  resistance  to  flow  between  capillaries  and  ventricular  CSF 
n  CF 

rCB=  resistance  to  flow  between  capillaries  and  brain  tissues 
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R 

CV  =  resistance  to  flow  between  capillaries  and  venous 
CCp  =  compliance  factor  between  choroid  plexus  and  ventricular  CSF 

£-(PC-Pf)  =  time  derivative  of  the  pressure  difference  between 
choroid  plexus  and  ventricular  CSF 

As  a  first  approximation ‘the  resistance  and  compliances  were 
considered  as  mean  effective  values  i.e.,  constants.  The  overall 
matrix  of  resistances  and  compliances  were  evaluated  via  an 
inverse  procedure  (Sorek  et  al .  1986a). 

Note  that  equation  (1)  describes  the  flux  balance  in  steady 
as  well  as  in  non  steady  situations.  The  cells  may  be  rigid 
giving  a  flux  term  expressed  by  pressure  differences  divided  by 
resistance  or  contractile  yielding  a  flux  which  is  the  product  of 
compliance  and  time  changes  of  pressure  differences. 

In  the  course  of  evaluating  the  resistances  of  the  model  it 
was  shown  that  a  situation  leading  to  evolution  of  a 
'normotensive'  hydrocephalus  may  take  place. 

The  usual  accepted  mechanism  for  development  of 
hydrocephalus  is  defective  absorption  of  CSF  in  the  venous  sinus. 

In  our  mathematical  model  this  will  be  expressed  by  RpS=°°  or 

^FS  =  1/R  =0  (where  R  __  is  the  resistance  to  flow  between  the 

r  b  FS 

ventricular  CSF  compartment  and  venous  sinus  compartment  and 

where  Zpg  stands  for  the  conductance  between  these  compartments). 

However  it  was  shown  (Sorek  et  al.  1986a)  that  at  the  same 
time  the  resistance  (R)  between  capillaries  (C)  and  brain  tissue 

(B)  attain  infinite  values  R„D  = 


o 
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In  the  set  of  equations  for  the  solution  of  the  model 
resistances  we  have  relations  between  fluxes,  pressure 
differences  and  coefficients  a  and  8  . 


a  = 


(2) 


8  = 


(3) 


For  instance  the  relation  for  (conductance 
capillaries  and  brain  tissue)  and  ZfS  (conductance 
ventricular  CSF  and  venous  sinus)  is  as  follows 


between 

between 


e(W 


CB  a<VPs)  •  (Pf"PB) 


(4) 


’”FS 


*  *  *  *  *  * 

a(pv-ps:>qF-(pF-yqA 

<yOi«<Vps>-<pF-V] 


(5) 


where  (  )*  denotes  mean  effective  values 

flux  entering  the  arterial  compartment 

n  flux  entering  the  CSF  compartment  =  CSF  generation 
WF 

QF  may  be  described  by 


*F 


P 


C 

R 


-P 


CF 


F 


(6) 


J 
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A  solution  of  the  set  can  be  attained  when  we  allow  e  =  0.  An 

accompanying  condition  Z  =0  will  still  yield  a  possible 

FS 

solution . 

Thus  the  following  mathematical  equations  now  exists 

ZCB=0  (Rcb  =  -  )  (7) 

zfs=o(rfs  =  »)  (8) 

an  #(P#  -P#  ).q  *(P#  _P  »)=0  (9) 

Equation  (7)  is  the  mathematical  representation  of  a  flow 
impediment  between  the  capillaries  and  the  brain  tissue.  The 

condition  Zcb_q  may  be  regarded  as  a  precondition  activating  the 
NPH  situation  represented  by  equations  (8)  and  (9). 

Conditions  expressed  by  equations  (8)  and  (9)  indicate 
blockage  of  CSF  transfer  from  the  ventricles  to  the  venous  sinus. 
As  this  conditions  do  not  affect  production  of  CSF  by  the  choroid 
plexus  (Qf=°)  compartment  (F)  will  upuna.  The  presence  of 
copliances  Cpg  c  and  C  (Fig.  3)  allow  for  the  expansion 

r  o  Lr 

without  increase  of  pressure  i.e.  NPH. 

Thus,  in  a  situation  where  the  flow  from  capillaries  to  brain 
tissue  is  impaired  as  may  be  the  case  arteriosclerotic 
cerebrovascular  disease  and  especially  in  small  vessel  disease  in 
the  aged,  a  NPH  may  develop. 

In  order  to  overcome  the  NPH  situation  one  has  to  interfere 
with  the  balance  as  stated  in  equation  (9).  By  lowering  the  CSF 

pressure  Pp  i.e.  shunting  procedure,  the  previous  conductivities 
may  change  in  a  step  fashion  accommodating  the  new  mean  pressures 


& 
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and  fluxes  as  indicated  by  equations  (4)  and  (5).  Note  that 
according  to  equation  (9)  such  step  changes  may  also  take  place 
-hen  changing  other  factors  e.g.  Pv,  Pji  atc.  Thus  removal  of 

CSP  will  also  yield  a  change  in  capillary  to  brain  tissue 
transfer  which  may  explain  the  improvement  in  neurological 
functions  after  shunting.  It  was  shown  that  CSF  drainage  in 
hydrocephalic  patients  increase  regional  cerebral  blood  flow 
(Symon  and  Hingzpeter,  1977).  If  the  small  vessel  disease 
continues,  equation  (9)  may  again  prevail  and  a  further  decrease 
in  CSF  pressure  is  necessary  in  order  to  accomodate  for  the  new 
resistances  and  complainces  as  could  be  the  case  in  our  patient. 
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FIG.  1  :  ENLARGED  LATERAL  VENTRICLES  INSPITE  OF  A  PATENT 
MEDIUM  PRESSURE  VENTRICULO  PERITONEAL  SHUNT. 


FIG.  2  : 


WELL  DRAINED  VENTRICLES  AFTER  INSTALLMENT 
OF  LOW  PRESSURE  V-P  SHUNT. 


=  750ml/min  Q,  =  750ml/min 


Lumped-parameter ,  eeven-compartmental  model  of  the  intracranial 
cerebrovascular  fluid  system  (R-resistance,  C-compliance ; 

C  3  pressure  (mm  Hg);  (  )  flow  (ml/min) ;  <  >  volume  (ml)). 
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A  SIMPLE  CONTINUUM  MODEL  OF  BRAIN  TISSUE  DEFORMATION 


by 


S..Sorek*,  J.  Bear**  and  Z.  Kami*** 


ABSTRACT 

The  continuum  approach  is  employed  for  modelling  the 
distribution  of  deformation  and  stress  within  the  brain  tissue 
visualized  as  the  domain  between  two  concentric  ellipsoids  having 
the  same  axes.  Phasic  pressures  in  the  Jugular  Bulb  and  CSF 
ventricles,  estimated  by  a  compar tmental  model,  provide  the 
conditions  on  the  outer  and  inner  boundaries,  respectively.  The 
simulation  yields  phasic  stress  and  deformation  which  are 
consistent  with  clinical  observations. 


Keywords:  brain  tissue,  CSF  ventricles,  jugular  bulb,  Hook’s 

law,  compartmental  phasic  pressures 
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INTRODUCTION 

Changes  in  pressure  in  both  the  jagular  bulb  and  in  CSF 
ventricles,  produce  stress  variations,  accompanied  by 
deformation,  within  the  brain  tissue.  Our  objective  in  this  paper 
is  to  obtain  estimates  of  stress  and  deformation  distribution 
within  the  brain.  To  achieve  this  goal,  a  simplified  continuum 
model  of  the  brain  is  introduced  and  investigated. 

In  the  literature  on  the  biomechanics  of  heart  impacts, 
much  use  of  spherical  and  oblate  spheroidal  coordinates  is  made. 
The  skull  is  idealized  as  a  rigid  sphere  with  an  opening  that 
stimulates  the  foramen  magnum,  while  the  spinal  dura  mater  is 
idealized  as  a  cylindrical  membrane  fitted  to  the  foramen  magnum 
(Liu,  1978).  The  intradural  content  of  the  central  nervous 
system  (CNS)  fluid-filled  continuum  is  regarded  by  some  as  a 
single-phase,  quasistatic  incompressible  fluid  in  a  two 
compartment  (skull  and  spinal  dura-mater)  structure,  or  as  a 
single-phase  compressible  elastic  fluid  possessing  a  single 
"averaged"  bulk  modulus  and  a  sheer  modulus  (Pollack  and  Boshes, 
1936;  Ommaya  and  Hirsch,  1971;  Lofgren  and  Zwetnow,  1973; 
Goldsmith,  1972;  Marmarou,  1973;  King  and  Chou,  1976).  Hakim  et 
al.  (1976)  stipulated  the  case  of  static  equilibrium  to  find  the 
stress  in  the  brain  tissue,  regarding  the  cranial  vault  as  a 
hollow  sphere  with  ventricles  in  the  middle.  The  inner  and  outer- 
boundary  conditions  were  the  pressure  value  of  CSF  in  the 
ventricles  and  at  the  subarachnoidal  space. 

Apart  from  the  spherical  models,  other  curvilinear 
mappings  were  suggested  such  as  a  truncated  shell  of  revolution 
(Schumacher,  1978). 
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Our  present  attempt  is  to  revive  Leonardo  da  Vinchi's 
idea,  based  on  his  anatomical  sketches  (Fig.  la)  (cf.  Russel, 
1959),  and  consider  the  domain  between  two  ellipsoids,  the  inner 
ventricular  ependyma  and  the  outer  pia  mater,  as  an  approximate 
geometry  of  the  .cranial  vault. 

The  brain  tissue  is  assumed  to  be  a  homogeneous  and 
isotropic,  single-phased  material. 

In  an  earlier  work,  Sorek  et  a  1  .  (1986a, b,c,),  simulated 
the  cerebrovascular  flow  regime  in  a  compartmental  model  cerebral 
system  (Fig.  2).  They  derived  the  non-steady  phasic  pressure  and 
flux  waves  associated  with  the  artery  (A),  capillary  and  choroid 
plexus  (C),  CSF  ventricles  (F),  veins  (V),  brain  tissue  (B), 
venous  sinus  (S),  and  jugular  bulb  (J).  Here  we  will  use  the 
phasic  pressures  in  the  CSF  ventricles  and  in  jugular  bulb 
obtained  in  that  study  as  an  inner  and  outer  boundary  conditions 
respectively,  imposed  on  the  considered  brain  tissue  domain. 

Clinical  data  show  that  the  material  comprising  of  this 
organ  is  deformable.  However,  it  is  not  strictly  elastic,  since 
under  external  stresses,  it  exhibits  also  the  delayed  property  of 
viscosity  (Pamidi,  1976).  Nevertheless,  as  a  first 
approximation,  we  choose  here  an  elastic  stress-strain 
constitutive  relation  for  this  material. 


104 


ASSUMPTIONS 

The  following  assumptions  underline  the  constructing  of 
the  model  simulating  stress  and  deformation  of  the  brain  tissue: 
1.  The  brain  tissue  configuration  is  approximated  as  a  shell 
between  two  prolated  spheroidal  bodies  (Fig.  lb).  Its  inner  (CSF 
venttricles)  and  outer  surfaces  are  obtained  by  rotating 
ellipses  about  their  common  major  axes.  The  relation  between  a 
cartesian  system  (X,Y,Z)  and  the  prolate  spheroidal  coordinate  (£., 
n,i|0  system  is  given  by 

x  =  a  sinh(£)  sin  (n)  cos  (i(i)  (1.1) 

y  =  a  sinh(C)  sin  (n)  sin  ( ip )  (l.<2) 

z  =  a  cosh(0  cos  (n)  (1.3) 

where,  ?>0,  0<n<tr;  (J<i(i<2ir  .  The  associated  scale  coefficients 

(h,  ,h  ,  h  ,)  are  expressed  by 
i  n  <• 

h2  =  h2  =  a2  [  sinh2(U  ♦  sin2(n)l  U.l) 

h2  =  a2  sinh2(C)  sin2  n  (t l  .d) 

The  equations  of  the  outer  and  inner  ellipsoids  are  given 
respectively  by  C  =  £„  =  const.  and  £  =  Cj  =  const . 


3.  The  brain  tissue  material  is  assumed  isotropic,  homogeneous 
and  elastic,  obeying  Hook's  Law,  which  can  be  expressed  by 


Jij  E  ^ij  °kk 


(3) 


where  w  -  Poisson's  ratio,  E  =  Young's  Modulus,  6i  j  =  Kronecker 
delta,  Eij  =  strain  tensor,  °ij  =  stress  tensor  and  0  i<k  =  d  i  agona  1 


stress  tensor. 
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4.  The  inner  and  outer  boundaries  are  loaded  by  spatially 
uniform  pressures  in  the  CSF  ventricles  and  in  the  jugular  bulb, 
respectively . 

5.  At  each  point  within  the  considered  domain,  we  can  determine 
the  radius  of  curvature  of  the  ellipsoid  C  =const.  passing 
through  that  point.  We  then  assume  that  deformations  are  mainly 
along  the  radius  of  curvature.  Also  it  is  assumed  that  the  organ 
is  undergoing  small  deformations.  The  compatible  relation  between 
strain  and  deformation  is  expressed  by 


-7  (u .  .  +u .  . ) 

2  1  ,j  j  ,1 


(4) 


where  u.  is  the  deformation  in  the  i  direction  and  u.  .  is  the 

1  i,J 

derivative  of  with  respect  to  the  j  coordinate. 

5.  Bending  moments  and  shear  forces  are  neglected  because  of  the 
symmetries  involved  in  the  assumed  configuration. 

6.  Body  forces  (e.g.,  due  to  gravity)  are  neglected. 


We  may  now  write  the  force  balance  equation  in  the  (4 , n * ^ ) 

coordinate  system.  Note  that  the  prevailing  stresses  are  normal 

to  the  surface,  o  ,  which  is  balanced  by  the  components  of  the 

longitudinal  3tress,  a  ,  and  the  latitudinal  stress 

nn 

FORCE  BALANCE  EQUATION 

Let  us  first  evaluate  the  length  of  the  mean  semimajor 

and  semiminor  axes  of  the  CSF  ventricles  visualized  as  a  prolate 

spheroidal.  Let  a^b^  and  cj  denote  the  semi-axes  of  the  Jagular 

bulb  and  a  ,  b  ,  c  denote  those  of  the  CSF  ventricles.  For 
F  F  F 

the  outer  surface,  we  introduce  the  estimates.  aj  =  9cm.;  bj=6cnt; 
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CJ ’ 

:3cm.  (  i  .  e . , 

a  volume 

of  approx 

i  m  a  t  e  1  y 

6  8  0  m  1.)  . 

The  volume  of  the  CSF 

ventricles 

is  taken 

as 

II 

u* 

> 

3 

150cm.  ,  i .e  . , 

V  =  -  la  be  = 

r  j  P  r  r 

150 

(5) 

We  now  assume  that  the 

inner  and  outer  prolate 

spheroids 

are 

similar.  This 

means  that 

bF  =  2CF 

(6) 

aF  ^  3CF 

(7) 

Hence,  by  virtue  of  equations  (5),  (6)  and  (7)  we  obtain 
a  =5  .4cm.  b  =3  .6cm .  c=1.8cm. 

r  r  r 

As  was  stated  above,  the  deformation  is  a  function  of 
the  K  coordinate  only.  Any  volume  element  (Fig.  3),  is 

subject  to  radial  Fr{.  ,  longitudinal  F  and  latitudinal  F 

nn 

forces.  The  lengths  of  the  volume  element  edges  are  h^,  hndn’ 

and  h  d  . 

*  * 

By  writing  a  force  balance  along  the  radius  of  curvature,  we 
obtain  the  resultant  in  that  direction  in  the  form. 

2F  sin  (4&  ♦  2F  sin  (4p)  *  dF  (8) 

♦♦  2  '  nn  2  U 


** 


For  snail  angles  d<i<<l  and  dn<<l  an  equivalent 

form  of  (8),  written  in  terms  of  stresses  is 
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—  (o..h,h  )  -  o,,h.h  -  o  h,h.  =  0 

d?  ec  <i»  n  ♦♦  £  n  nn  £  (g) 

One  may  obtain  (9)  by  replacing  the  forces  in  (8)  by  stresses 
multiplied  by. the  appropriate  areas. 

Equation  (9)  is  subject  to  the  following  conditions  on  the 

external  (4  =  E  )  and  internal  (  £  =  £  )  surfaces. 

J  F 


O  s  _P 

££  pf 


at  £  =  £ 


(10.1) 


°££  -  -PJ 


£  =  £ . 


(10  .2) 


We  note  that  because  of  symmetry,  the  resultants  in  the 
other  two  directions  vanish  . 

Recalling  the  assumption  that  the  displacement  u,  is  normal 
to  the  £  ^constant  surface  (i.e.,  u  =  u(£))  and  that  these 
displacements  are  very  small,  we  now  stipulate  the  condition 


u=o  at  5  =o 


In  view  of  equations  (1 ),(2),  (4)  and  (11),  we  obtain 


'££  d£ 


e  ■ 

nn 


o  Ihn(£+u)~  hn(£)ldn  _  (H-du/d£)2(S+u)2-£2 


/  hn(ndn 


-  2  -  • 
*  2  £  ’ 


«  l) 


(13) 
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I  lh»(C*-)'b*<5>ld*  _  1>du/d6)(£,ul.£  u 

(14) 

it 

l  w 

?  '  e 

We  Mill  relate  strain  to 

(equation  (3))*  By  virtue 

therefore  obtain 

stress  by  Hook's 

of  equations  (3) 

constitutive 

and  (l2)-(Hl) 

dc  "  E  l°ec_M(0nn+a4'V)  J 

( 15 ) 

2!’f  '“nn'1*1  V°£t)l 

(16) 

—  =  —  I  (j  —  u  ( a  +o  )} 

C  E  1  <14  ^  K  nn  J 

(17) 

Equating  equations  (16)  and 

a  =  Ao  +  Bo 

<|n|»  nn  ££ 

where 

(17),  we  obtain 

(18) 

A  =i±2E 

2+n 

(19.1) 

2+w 


( 19  .2) 


Substituting  equation  (18)  into  (9)  yields, 


nn 


d0« 

c  -~r  *  Do 

a  t 


(20) 


•there 


c 


VAh 


n 


(21 .1) 
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D 


■jr  (h  h  )-Bh.h 
dg  4»  n  Co 

h.(h ,+Ah  ) 

£  v  n 


(21 .2) 


Differentiating  (16)  with  respect  to  £  ,  equating  the  result  to 
(15),  substituting  (18)  and  (20),  we  obtain 


(22) 


where 


M  =  c  C 


(23.1) 


(23.2) 

(23.3) 


We  then  solve  equation  (22)  together  with  boundary 
conditions  (10),  for  =  o^(£)  .  With  the  calculated  o ^  , 

we  evaluate  ,  o  and  u  as  functions  of  £  ,  using 

equations  (20),  (18)  and  (17)  respectively.  All  calculations  are 

done  numerically. 
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IMPLEMENTATION 

Employing  equations  (1),  with  the  values  a,b,  and  c  of 
the  semi-axis,  of  the  inner  and  outer  surfaces,  we  calculate 
5  =  0.454,  c j=  0.689  and  ^=^=26.5°,  atnF3ij=90o.  We  then 

impose  the  phasic  pressures  Pp  and  Fj  ,  evaluated  by  Sorek  et  al. 
(1986),  using  compartmental  modelling,  as  boundary  conditions 
(equations  (10))  to  calculate  phasic  variations  of  stressess  and 
deformation  as  distributed  along 

Figure  4  depicts  the  mean  distribution  of  stressess  arr  , 0 

tin 

,0^  and  deformation  u  due  to  mean  pressures  PF=10mmHg  and 
Pj=2mmHg.  Note  that  compression  stress  in  result  in  tension 

in  the  longitudinal  ,  0  ,  and  latitudinal  ,  ,  stresses. 

Also  note  the  hyperbolic  characteristic  for  deformation  and 
stresses,  demonstrating  a  decay  in  intensity  from  the  CSF 
ventricles  surface  to  the  outer  brain  tissue  surface.  Figures  5 
and  6  show  0^  ,  and  u  surfaces  as  functions  of  time  and 
location,  depicting  the  pulsational  nature  of  stresses  and 
deformation.  The  range  of  deformation  as  described  by  figure  6 
is  consistent  with  reported  clinical  observation  (Allen  et  al. 
1983) . 

CONCLUSION 

A  continuum  model  was  developed  and  employed  to  simulate 
stress  and  deformation  in  an  elastic  ellipsoidal  shaped  brain 
tissue.  The  results  of  a  compartmental  model  were  employed  to 
account  for  phasic  pressures  used  as  boundary  conditions  in  the 
model  presented  herewith.  Although  the  model  is  basically  a  1-0 
and  simplified  results,  are  consistent  with  clinical 


observations.  Modifications  such  as  the  use  of  a  non-elastic 
stress-strain  constitutive  la-,  may  yield  better  predictions  of 
stress  and  deformation  distribution  in  time  and  space. 
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